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Abstract 


This thesis presents a study of the effect and generation of non-adiabatic pertur¬ 
bations in Cosmology. 

We study adiabatic (curvature) and entropy (isocurvature) perturbations pro¬ 
duced during a period of cosmological inflation that is driven by multiple scalar 
fields with an arbitrary interaction potential. A local rotation in field space is 
performed to separate out the adiabatic and entropy modes. The resulting field 
equations show explicitly how on large scales entropy perturbations can source 
adiabatic perturbations if the background solution follows a curved trajectory in 
field space, and how adiabatic perturbations cannot source entropy perturbations 
in the long-wavelength limit. It is the effective mass of the entropy field that de¬ 
termines the amplitude of entropy perturbations during inflation. We show why 
one in general expects the adiabatic and entropy perturbations to be correlated 
at the end of inflation, and calculate the cross-correlation in the context of a 
double inflation model with two non-interacting fields Jl|. Then, we consider 
two-field preheating after inflation, examining conditions under which entropy 
perturbations can alter the large-scale curvature perturbation and showing how 
our new formalism has advantages in numerical stability when the background 
solution follows a non-trivial trajectory in field space [||, ||. 

Then we compare the latest cosmic microwave background data with theoret¬ 
ical predictions including correlated adiabatic and CDM isocurvature perturba¬ 
tions with a simple power-law dependence. We find that there is a degeneracy 
between the amplitude of correlated isocurvature perturbations and the spectral 
tilt. A negative (red) tilt is found to be compatible with a larger isocurvature 
contribution. The main result is that current microwave background data do not 
exclude a dominant contribution from CDM isocurvature fluctuations on large 
scales, and marginally favour a significant fraction |Q . 

We then study perturbations in Randall-Sundrum-type brane-world cosmolo¬ 
gies. The density perturbations generate Weyl curvature in the bulk, which in 
turn backreacts on the brane via stress-energy perturbations. On large scales, 
the perturbation equations contain a closed system on the brane, which may be 
solved without solving for the bulk perturbations. Bulk effects produce a non- 
adiabatic mode, even when the matter perturbations are adiabatic, and alter the 
background dynamics. As a consequence, the standard evolution of large-scale 
fluctuations in general relativity is modified. The metric perturbation on large- 
scales is not constant during high-energy inflation . 

The effect of non-linear perturbations on initiating inflation is examined from 
the perspective of both spacetime embedding and scalar field dynamics. Scalar 
field dynamics that is consistent with the embedding constraints are examined, 
with the additional treatment of damping effects. The effects of inhomogeneities 


on the embedding problem also are considered. A category of initial conditions 
are identified that are not acausal and can develop into an inflationary regime 

i- 

Finally the work is summarised and current and future extensions are dis¬ 
cussed. 
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Chapter 1 
Introduction 


In this thesis I have aimed to present the work I did during my PhD. It is based on 
the publications and a preprint I collaborated on during that time [|l], . 

My other publication |P] was completed before starting my PhD and so is not 
included in the thesis. 

In this introduction Chapter I try to place the work done during the thesis 
in context. To do this I give a summary of the current state of cosmology both 
in terms of theory and experiment. My summary is centered around the issue of 
adiabatic and entropy perturbations in Cosmology and is by no means exhaustive. 
There are a large number of good text books available where a more exhaustive 
treatment is available, see for example []^ 


1.1 Big bang cosmology 


The big bang model is based on the assumption that the Universe is isotropic 
and homogeneous on large scales. There is much evidence for this. In particular 
the cosmic microwave background (CMB) has been measured to be isotropic to 
one part in 10® |^. Together with the weak Copernican principle, i.e. all cosmic 
observers see a nearly isotropic CMB, this implies that the universe is nearly 
homogeneous on large scales [M. Large scale structure studies support this. 


This assumption then leads to the following space-time metric: 


= —dt^ -I- a{ty 


dr'^ 


1-Kr^ 


-h sin^ 9 + d(p^ 


where a{t) is the scale factor, t is the time coordinate, r, 6, (p are the spatial polar 
coordinants and K is the curvature which can be negative, zero or positive. 

Using the Einstein equations and assuming the matter in the universe is a 
mixture of perfect fluids we get the Friedman equation: 


Hitf 


SttG , , K 


( 1 . 1 ) 
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where H = a/a is the Hubble parameter, G is Newton’s constant, p is the energy 
density. 

In the 1920’s Hubble observed that the red shift of light emitted from galaxies 
increases with their distance and thus showed that H is positive and so the 
Universe is expanding. If a is positive then a must have been zero at some time 
in the past. This model of a Universe expanding from an initial singularity is 
often referred to as the ‘big bang’ model. However at very high energy densities 
quantum corrections are thought to change the Friedman equation and perhaps 
avoid the singularity. 

From energy momentum conservation we get the continuity equation 

d{pa^) = —pd{a^) 

where p is the pressure. For the simple equation of state 


p = wp, 


we then get 


p (X a 


3(l+-!i>) 


Some cases of interest are 


Radiation: 

(P = \p) - 

-s- p oc 

Matter: 

(p = 0) - 

-s- p oc a“^. 

Vacuum energy: 

(p = -p) - 

-s- p oc constant 


So for a Universe consisting of matter, radiation and vacuum energy (also known 
as a cosmological constant) the Friedman equation becomes 

2 SttG fpr{to)a{toY pm{to)a{tof , K 

=~( aW a{€f 

where subscripts r, m and v refer to radiation, matter and vacuum respectively. 
The constants of integration are set at some initial time to- As a is increasing 
with time, it can be seen from the above equation that the Universe can go 
through different stages where different components in the Friedman equation 
will dominate. Clearly if pv{to) is non-zero and K is non-positive, then eventually 
the vacuum energy will dominate no matter how large the other initial densities 
or curvature are. If K is positive then it is possible that the expansion will be 
halted and reversed before the vacuum energy dominates the expansion. Current 
experimental data roughly indicates the following sequence of events: About 
15 billion years ago there was the big bang, and the Universe was radiation 
dominated. After several hundred thousand years, the Universe became matter 
dominated. Then about three billion years ago (red-shift one), the Universe 
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became vacuum energy dominated. In order for the Universe to have such a long 
period of radiation and matter domination, K must be very close to zero. 

For radiation the temperature T is related to the energy density by 

p oc T^. 

Today the radiation component of the Universe is still detectable even though it 
only contributes a minute fraction of the current energy density of the Universe. 
It has been red-shifted to microwave frequencies and is known as the cosmic 
microwave background. It has a black body spectrum due to the radiation being 
in thermal equilibrium. Its detection was one of the main pieces of evidence for 
the big bang theory. Thus the temperature magnitude and evolution of the early 
Universe can be inferred. This allows the modelling of the formation of helium and 
the other light elements from hydrogen during the early Universe, a process known 
as “big bang nuclear synthesis” (BBN). These inferred proportions of hydrogen, 
helium and other light elements are in good agreement with observations and 
provide further solid evidence for big bang Cosmology. 

1.2 Structure formation 

Although the Friedman equation provides an excellent description of the be¬ 
haviour on average, there is still the task of explaining the formation of galax¬ 
ies and galaxy clusters. Clearly gravitational attraction will cause any inhomo¬ 
geneities present in the Universe to increase with time. So the early Universe 
would have to be more homogeneous than today, but still there need to be some 
inhomogeneities that can grow into the large scale structure we see today. 

In the hot early Universe, the CMB radiation would have been coupled to the 
baryons by Thomson scattering. The average path lengths of the photons would 
have been very short as a consequence. However as the Universe expanded and 
the temperature dropped to about 10^ K the radiation became decoupled from 
the baryons and so the path length of the photons became almost unlimited. 
Thus, since several hundred thousand years after the big bang, the radiation 
photons would have been free streaming. As a consequence when we observe 
them, we are able to see a snap shot of what the Universe looked like at the time 
of decoupling. In effect we have a picture of the spatial temperature variation on 
a sphere surrounding us with a radius of about hfteen billion light years. This 
allows us to infer what the density inhomogeneities were at that time. 

The hrst accurate detection of the inhomogeneities (or perturbations) was 
made by the COBE satellite [Q . It measured the temperature variation on degree 
scales and showed it to be one part in 10^. It also found that the power spectrum 
of the perturbations was close to scale invariant. This is a particularly simple 
perturbation spectrum as it implies there is no characteristic length scale for the 
perturbations. 
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In reality the scale invariance is broken on scales smaller than COBE was 
able to measnre. This is because the equations describing the evolution of cos¬ 
mological perturbations (see for example [11]) show that a characteristic scale 


for the perturbations is the Hubble length 1/H. Typically the adiabatic pressure 
does not effect the perturbations on scales larger than the Hubble length. During 
the standard cosmological evolution the Hubble radius is increasing in size faster 
than the scale of the perturbations, i.e. for a perturbation with a wave number 
k/a, the ratio k/{aH) will be increasing with time. Thus perturbations with 
scales initially greater than the Hubble length will eventually have scales smaller 
than the Hubble length. At this stage, pressure can counteract the tendency of 
gravitational collapse of the perturbations. This can lead to ‘acoustic oscillations’ 
where the magnitude of the perturbation oscillates with time. These oscillations 
lead to peaks in the CMB power spectrum which have been observed in recent 
CMB experiments [O, O, O. 


1.3 Inflation 

Although the big bang model is well validated by experiment, it does seem to 
require rather special initial conditions. At the time of decoupling the Hubble 
horizon (which is roughly equal to the causal horizon) would only be about one 
square degree. Yet the CMB is the same temperature to one part in 10® in all 
parts of the sky. What’s more the deviations from homogeneity have the special 
form of having a scale invariant power spectrum on scales larger than the Hubble 
horizon. In the big bang model there is no causal way of achieving this special 
setup and it just has to be put in by hand as an initial condition. These are known 
as the ‘homogeneity’ and ‘inhomogeneity’ problems, i.e. why is the Universe so 
homogeneous and why are the deviations from homogeneity of such a special 
form. 

Another unusual thing is how close the curvature K must be to zero in order 
to have such a long radiation and matter dominated era. This is known as the 
flatness problem. 

A natural solution to these and various other problems can be provided by 
the model known as inflation (see for a recent review). From particle physics 
models one expects that under the extreme conditions of the early Universe there 
would be one or more scalar helds. It can be shown (see for example [|^) that pro¬ 
vided the potential energy of the scalar held dominates the gradient and velocity 
energy in a volume larger than a Hubble volume and the potential is sufficiently 
hat, then any inhomogeneities in the held will be smoothed out. The criterion 
that the potential energy has to dominate over a Hubble volume to some extent 
undermines the solution to the horizon problem. If the ‘chaotic inhation’ sce¬ 
nario is considered, where the scalar held is taken to be distributed randomly 
throughout the Universe at the Planck time, then it can be shown that in some 
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areas this criterion will be met just by chance, i.e. there is a hnite probability 
of the criterion being met and so if there is a large area of randomly fluctuating 
scalar held, there is a large probability that the criterion will be met somewhere. 
However if one wants to consider inhation starting in some isolated patch well 
after the Planck time, such as in ‘new inhation’, then the homogeneity criterion 
becomes more restrictive. This issue is examined further in Chapter p. 

The other criterion that the potential must be sufficiently hat is quantihed by 
the slow roll conditions: 



( 1 . 2 ) 


(1.3) 


where 0 is the scalar held value, primes denote diherentiation with respect to 0 
and V denotes the potential of 0. It can be shown that if the slow roll conditions 
are satished then the scale factor accelerates with time and the curvature is driven 
to zero and so inhation solves the hatness problem. To be more specihc 


a > 0 


where dot denotes diherentiation with respect to proper time. This is in contrast 
to ordinary matter where a deaccelerates and the curvature moves away from 
zero. 

For a homogeneous scalar held, the evolution is given by the Klein Gordon 
equation 


0 + 3iP0 H—j— — 0 
dcp 


For a homogeneous scalar held dominated Universe the Friedman equation is 
given by 



Applying the slow roll conditions, the Klein Gordon and Friedman equations can 
be approximated by 



and 


3if0^ -U'(0). 
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A common potential used is 


V = 

2 

where \/V^ = m is the effective mass of 0 . 

Remarkably inflation also solves the inhomogeneity problem. On sub-Hubble 
scales, quantum fluctuations in the value of 0 occur. In inflation the Hubble 
scale is almost constant, while the scale factor is rapidly increasing. Thus the 
sub-Hubble vacuum fluctuations eventually are stretched to super-Hubble scales. 
At this stage they lead to classical perturbations in the value of 0 on super-Hubble 
scales. 

It can be shown (see for example |^) that the power spectrum of fluctuations 
generated by inflation has spectral index 

n = 1 + 27] — 6e 

where n = 1 represents scale invariance. These fluctuations are inherited by 
the radiation era when the scalar field decays into radiation and other forms of 
matter after inflation ends. So remarkably, the same slow roll conditions solve 
the homogeneity, flatness and inhomogeneity problems. 

1.4 Reheating 

Inflation ends when the slow roll conditions are violated. This is the stage when 
0 (the ‘inffaton’) oscillates in the valley of the potential V. It is then thought to 
decay into photons ( 7 ), non-baryonic or ‘cold dark matter’ matter (c), baryonic 
matter ( 6 ) and neutrinos (z/). This process is known as ‘reheating’. 

Simple aspects of this process of the transfer of energy from the inffaton can 
be modelled by the addition of a scalar field y which has a coupling term 
in the potential, where g is a. coupling constant. Neglecting perturbations in the 
metric (a full discussion or perturbations is given in Chapter the equation of 
motion for the perturbation of wave number A; of y is given by: 

Sxk + 3H6xk + SXk = 0 . 

At the end of inflation 0 is oscillating in the valley of the potential. For V = 
the behaviour of 0 is well approximated by 

0 ~ 00 sin mt 

where 0 o is the amplitude of the oscillations, for small wave lengths {k/a S> 
H) the expansion can be neglected. The perturbation equation can then be 
transformed to the well known Mathieu equation (see for example ||16|| ) 

+ (^ - 2,qcos2z)5xk = 0 
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where A = /{m^o?) + 2g, q = 5f^0g/(4m^), and z = mt. The behaviour of this 

equation can be understood with the aid of the Mathieu chart where there 
exist bands of instability where for certain value of A and q there is exponential 
growth of 5xk- Including expansion means that modes of 5xk move in and out of 
bands of resonance. This process is known as ‘preheating’ 0 . 0 - This represents 
the transfer of energy from the homogeneous held (j) into the perturbations of x- 
The homogeneous held (j) represents a condensate of zero momentum particles 
while 6xk represent bosons with momentum k. Thus preheating models the 
exponential production of bosons. Fermions can also be considered. Once most 
of the energy of 0 has been transferred into 6x then the y and remaining 0 
particles can in turn decay into other particles which eventually thermalize and 
make up the contents of the radiation era. 

The ehect of having a preheating phase is to end the oscillating phase of 0 more 
rapidly and it can lead to a higher temperature. There is also the possibility that 
the super-Hubble scale perturbations could be amplihed. This would have the 
undesirable consequence of the complicated and uncertain physics of preheating 
having an effect on the observed anisotropies in the CMB. This issue will be 
discussed further in Chapter 


1.5 Adiabatic and entropy perturbations 

Another special property of the fluctuations generated by a single scalar held is 
that they are adiabatic, i.e. the perturbation in 0 can be expressed as a time 
shift in the background scalar held: 

6(j) = ^{t)6t{t, x) 

where x is the spatial position vector and 6t is the time shift of the background 
solution needed to produce the perturbation. 

The adiabatic nature of the perturbations is inherited by the perturbations 
in the radiation era once the scalar held decays. So we have 

_ ^pu ^pc ^ pb 

P'y Pu pc Pb 

Using the continuity equation 

p = -3H{p + p) 

the adiabatic condition becomes 

3 ^ P'y 3 S Pi/ ^ Pc ^ Pb 

P'y Pu Pc Pb ' 
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This can also be expressed in terms of perturbations in particle number n 
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= 0 


where i and j can be 7, z/, c or b. Although the adiabatic condition is quite re¬ 
strictive, it does seem to be in good agreement with observations (see for example 
) and as mentioned above is also predicted by the simplest inflation¬ 
ary models. However a priori there seems to be no reason why the primordial 
perturbations have to satisfy the adiabatic condition. 

For example if there are two scalar helds (0 and x) present during inflation 
then the perturbations in each held do not necessarily have to correspond to the 
same time shift, i.e. 


12, 13, 14 


0 ^ X ' 

The generation of and evolution of adiabatic perturbations during multi-held 
inhation models is examined further in Chapter ^ Another important issue that 
is discussed is the possibility of there being correlations between the adiabatic 
and entropy perturbations. 

If there are entropy perturbations in the inhation stage, these can lead to 
entropy perturbations in the radiation dominated phase. For example we could 
have 

bpc ^ 3 b 
Pc ^ 4: P^' 

The entropy perturbation is dehned as that part of the perturbation which does 
not satisfy the adiabatic condition, i.e. 

(o ^ Pc 3 6p-f 

‘^C'Y T 

pc 4 

where Sc-y is the entropy perturbation between cold dark matter and photons. 

The ehect of entropy perturbations in the radiation era on the observed CMB 
spectrum is discussed in Chapter Using current data the magnitude and cor¬ 
relation of the entropy perturbation is estimated. 

1.6 Brane world cosmology 

String/M-theory theory predicts the existence of extra spatial dimensions. As 
a simple example, our Universe could be a three plus one hyper-surface (the 
‘brane’) embedded in a four plus one space-time (the ‘bulk’). 
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Until recently it was thought that any extra dimensions would be too small to 
be detectable, since it was thought that the deviations in gravity would contradict 
existing experiments. However Arkani-Hamed et. al. [113 presented models with 


large compact extra dimensions, and then Randall and Sundrum (RS) ^ 


showed that even an inhnite extra dimension is possible since gravity can be 
localized near the brane by the curvature of the bulk. 

In the RS model, the bulk is anti-de Sitter and the brane has positive vacuum 
energy (the ‘tension’) which cancels out the bulk negative vacuum energy. The 
particles of the standard model are conhned to the brane. 

There has been much interest in seeing whether such a scenario would have 
any cosmological implications (see for a review). In particular the question 


of whether the extra dimension would lead to an identihable signature on the 
CMB is of great interest. 

The bulk affects the perturbation equations on the brane in two ways. It adds 
terms which are quadratic in the brane energy momentum tensor. It also adds 
terms from the projected Weyl tensor of the bulk. These can be viewed as an 
additional fluid on the brane and so lead to the possibility of additional entropy 
perturbations. This and other related issues are discussed in Chapter |. 
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Chapter 2 

Adiabatic and entropy 
perturbations from inflation 


2.1 Introduction 


As discussed in the introduction, inflation in the early universe has become the 
standard model for the origin of structure. Inhomogeneities in the present matter 
distribution can be traced back to quantum fluctuations in the fields driving 
inflation which are stretched beyond the Hubble scale during inflation. In the 
simplest models of inflation driven by a single scalar field, these fluctuations 
produce a primordial adiabatic spectrum whose amplitude can be characterized 
by the comoving curvature perturbation TZ, which remains constant on super- 
Hubble scales until the perturbation comes back within the Hubble scale long 
after inflation has ended. 

As soon as one considers more than one scalar held, one must also consider 
the role of non-adiabatic fluctuations. This can have important consequences, 
both in affecting the evolution of the curvature perturbation (often referred to as 
the ‘adiabatic perturbation’), but also in the possibility of seeding isocurvature 
(or ‘entropy’) perturbations after inflation. 

Previous studies have demonstrated that non-adiabatic pressure perturbations 
can alter the curvature perturbation on super-Hubble scales either during infla- 
or after 1^ 


tion [^, 


33. A general formalism to evaluate the 


curvature perturbation at the end of inflation in multiple field models was devel¬ 
oped in Ref. |^. In the presence of non-adiabatic fluctuations, one must follow 


the evolution of perturbed fields on super-Hubble scales, in particular tracking 
the perturbation in the integrated expansion [^, ^ 


in order to 


evaluate the large-scale curvature perturbation at late times [§§ 0 0, 


33, 34, 39 


However no similar formalism has been developed so far to evaluate the isocur¬ 
vature perturbation in the general case. Instead, isocurvature perturbations have 
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been studied in a number of particular models of inflation . These fluc¬ 

tuations typically arise as baryon modes (e.g. |^) or cold dark matter modes 
but neutrino isocurvature modes have also been considered [Q. Recently, it has 
been pointed out ^ that it is rather natural to expect the curvature and 


isocurvature perturbations to be correlated, which yields distinctive observational 
results [^, in contrast to the pure or uncorrelated isocurvature perturbations 
usually tested against observations |4F 


In this Chapter we will develop a general formalism to study the evolution 
of both curvature and isocurvature perturbations in a wide class of multi-field 
inflation models by decomposing field perturbations into perturbations along the 
background trajectory in field space (the adiabatic field perturbation), and or¬ 
thogonal to the background trajectory (the entropy field). We allow an arbitrary 
interaction potential for the fields, and, although we concentrate upon the case 
of two scalar fields, the general approach can be easily extended to N fields, 
where there will be iV — 1 entropy fields orthogonal to the background trajectory. 
This was done for a specific assisted inflation model in Ref. [TO. We will work 


in the metric based approach of Bardeen in order to define gauge-invariant 
cosmological perturbations, but our formalism can also be applied to the study 
of multiple scalar fields in other approaches |^, |^, . 


We begin by reviewing the standard results obtained in single field models, 
emphasizing the suppression of non-adiabatic fluctuations on large-scales. We 
then extend our analysis to general two-field models, defining an adiabatic field 
and an entropy field, whose fluctuations, though uncorrelated on small scales, 
may develop correlations through the subsequent evolution. 


2.2 Perturbation equations for multiple scalar 
fields 

We consider N scalar fields with Lagrangian density: 

1 ^ 

c = -V (Lfi, ■ ■ ■ , lpn) - , (2.1) 

^ i=i 

and minimal coupling to gravity. In order to study the evolution of linear 
perturbations in the scalar fields, we make the standard splitting ipj{t,x.) — >• 
(flit) + Sipi{t,x.). The field equations, derived from Eq. ( p.l|) for the background 
homogeneous fields, are 


(pi + SHifj Vtpj — 0 , 


( 2 . 2 ) 
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where = dV/dx^ and the Hubble rate, if, in a spatially flat Friedmann- 
Robertson-Walker (FRW) universe, is determined by the Friedman equation: 


= 



(2.3) 


with a{t) the FRW scale factor. 

Consistent study of the linear held fluctuations Sipi requires that we also con¬ 
sider linear scalar perturbations of the metric, corresponding to the line element^] 

ds^ = —{1+ 2A)dt^+ 2aB^idx''dt 

+ a? [(1 — 2%jj)5ij + 2E^ij] dx^dx^ , (2.4) 


where we have not at this stage specified any particular choice of gauge |^, ^ Q. 

Scalar held perturbations, with comoving wavenumber k = 2na/X for a mode 
with physical wavelength A, then obey the perturbation equations 

6(pi 3H6(pj H— 2^^^ ^ ^ VipjipjSifij 


— —2V^jA -|- if I 


■ k"^ ■ 

H -I- 3-^ H—— aB) 


(2.5) 


The metric terms on the right-hand-side, induced by the scalar held perturba¬ 
tions, obey the energy and momentum constraints 


..2 r 


3H [jj + HAj + ^ [^ + H{a^E - aB) 

ip +HA 


—AirGSp , 
—4:7rG6q. 


( 2 . 6 ) 

(2.7) 


The total energy and momentum perturbations are given in terms of the scalar 
held perturbations by 


6p = [f/ 

1 

dq,i = -Y ■ 

I 


( 2 . 8 ) 

(2.9) 


These two equations can be combined to construct a gauge-invariant quantity, 
the comoving density perturbation 


Cm = dp - 3H6q 


I 


\pi - viA - 'ipiS'-pi 


( 2 . 10 ) 


^ We follow the notation of Ref. |11 , apart from our use of A rather than (p as the pertur¬ 
bation in the lapse function. 
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which is sometimes used to represent the total matter perturbation. 

Because the anisotropic stress vanishes to linear order for scalar helds mini¬ 
mally coupled to gravity, we have a further constraint on the metric perturbations: 


a^E 


aB] +H 


a^E 


aB) — ^4 = 0 . 


( 2 . 11 ) 


The coupled perturbation equations (|2.5|) - (|2.9|) and (|2.11|) are probably most 
often solved in the zero-shear (or longitudinal or conformal Newtonian) gauge, 
in which a'^Ef — aBg = 0 |^. The two remaining metric perturbation variables 


which appear in the scalar held perturbation equation, = $ and are 

then equal in the absence of any anisotropic stress by Eq. (|2.11|). 

a 


Another useful choice is the spatially hat gauge, in which ^pQ = 0 


The 

scalar held perturbations in this gauge are sometimes referred to as the Sasaki or 
Mukhanov variables |Q, which have the gauge-invariant dehnition 


Qi = • 


( 2 . 12 ) 


The shear perturbation in the spatially hat gauge is simply related the curvature 
perturbation, T, in the zero-shear gauge: 

a'^Eq - aBq = a^E - aB + . (2.13) 

H H 

The energy and momentum constraints, Eqs. (|2.6| ) and (|2.7|) , in the spatially hat 
gauge thus yield 


— \ 1 > = -iTiGer. 
a^ 


(2.14) 

HAq = -iirGSqQ, (2.15) 

where is given in Eq. (p.lO|) , and from Eq. ( |2.9| ) we have Sqq = — ViQi- 


The equations of motion, Eq. ( p.5|) , rewritten in terms of the Sasaki-Mukhanov 
variables, and using Eqs. ( 2.14 ) and ( |2.15|) to eliminate the metric perturbation 
terms in the spatially hat gauge, become 


Qi + 3HQj H —-Qj 




V. 


SttG 






Qj = o. 


(2.16) 


2.2.1 Curvature and entropy perturbations 


The comoving curvature perturbation [^, ^ is given by 




p) - 


H 


P + P 


5q 


E 


Pi 


Qi ■ 


(2.17) 
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This can also be given in terms of the metric perturbations in the longitudinal 
gauge as 


11 


H 

7^=T-^(T + T^<^> 
H 


(2.18) 


For comparison we give the curvature perturbation on uniform-density hyper¬ 
surfaces, 


-C,=^ + H 


5p_ 

P 


(2.19) 


hrst introduced by Bardeen, Steinhardt and Turner as a conserved quantity 
for adiabatic perturbations on large scales ^1], It is related to the comoving 
curvature perturbation in Eq. (|2.17D by a gauge transformation 


-C = 7^ + 


2p 


3(p + p) \aH 


k 




( 2 . 20 ) 


where we have used the constraint equation ( p.l4| ) to eliminate the comoving 
density perturbation, Cm- Note that 71 and —( thus coincide in the limit /c —> 0. 

Both 71 and —( are commonly used to characterise the amplitude of adia¬ 
batic perturbations as both remain constant for purely adiabatic perturbations 
on sufficiently large scales as a direct consequence of local energy-momentum con¬ 
servation , allowing one to relate the perturbation spectrum on large scales to 
quantities at the Hubble scale crossing during inflation in the simplest inflation 
models |62| . 

A dimensionless definition of the total entropy perturbation (which is auto¬ 
matically gauge-invariant) is given by 


S = H 


6p 6p 

P P 


( 2 . 21 ) 


which can be extended to dehne a generalised entropy perturbation between any 
two matter quantities x and y: 


^xy I . 


X 


( 2 . 22 ) 


The total entropy perturbation in Eq. ( |2.21|) for N scalar helds is given by 
2 (h + E j - PI A) 


S = 


3(2I/ + 3ifE,7 0^j) 


(2.23) 


where the perturbation in the total potential energy is given by 5V = 'Yli 


14 




























The change in TZ on large scales (i.e., neglecting spatial gradient terms) can be 


directly related to the non-adiabatic part of the pressnre pertnrbation |2^, ^ 


n ^ -3H-S. 
P 


(2.24) 


We will now consider the evolntion of the adiabatic and entropy pertnrbations in 
both one- and two-field models of inflation. 


2.2.2 Single field 

Pertnrbations in a single self-interacting scalar field obey the gange-dependent 
eqnation of motion 


(2.25) 


6ip 3H6ip -\- I — -|- Vip^p 1 Sip 


— —2V^A -p ip 


■ ■ 

A + 3'ip — aB) 


snbject to the energy and momentnm constraint eqnations given in Eqs. 

The scalar field pertnrbation in the spatially flat gange has the gange-invariant 
definition, Eq. ( p.l2| ), 


ip 

Qif = Sip + —"0 . 


(2.26) 


For a single field this is directly related to the cnrvatnre pertnrbation in the 
comoving gange, where the momentnm, Sq = —ipSip, vanishes 


H H 

K = 'il) + —Sip = —Q^ 

ip ip 


(2.27) 


It is not obvious that the intrinsic entropy perturbation for a single scalar 
field, obtained from Eq. ( p.23|) . 


5 = 


2W 


3ip‘^{3Hip + 2VA) 


ip [Sip — ipA ) — ipSip 


(2.28) 


should vanish on large scales. Because the scalar field obeys a second-order 
equation of motion, its general solution contains two arbitrary constants of inte¬ 
gration, which can describe both adiabatic and entropy perturbations. However 
S for a single scalar field is proportional to the comoving density perturbation 
given in Eq. (|2.10|) , and this in turn is related to the metric perturbation, T, via 


Eq. (|2.14|) , so that 


5 = 




67rGipA3Hip + 2V^ 


—T 
a2 


(2.29) 
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In the absence of anisotropic stresses, \E' must be of order Aq, by Eq. (|2.11|) , 
and hence the non-adiabatic pressure becomes small on large scales [0, 

The amplitude of the asymptotic solution for the scalar held at late times (and 
hence large scales) during inhation thus determines the amplitude of an adiabatic 
perturbation. 

The change in the comoving curvature perturbation is given by 

(2.30) 

and hence the rate of change of the curvature perturbation, given by d In T^/d In a ~ 
(fc/aid)^, becomes negligible on large scales during single-held inhation. 

2.2.3 Two fields 

In this section we will consider two interacting scalar helds, cf) = ipi and y = 
if 2 - The analysis developed here should be straightforward to extend to include 
additional scalar helds, but we do not expect to see any qualitatively new features 
in this case, so for clarity we restrict our discussion here to two helds. 

In order to clarify the role of adiabatic and entropy perturbations, their evo¬ 
lution and their inter-relation, we dehne new adiabatic and entropy helds by a 
rotation in held space. The “adiabatic held”, cr, represents the path length along 
the classical trajectory, such that 

d = (cos6')(;/)-|-(sin^jy , (2-31) 

where 

cos^ = — , sin6'=— ^ (2.32) 

This dehnition, plus the original equations of motion for 0 and y, give 

a + 3Ha + I4 = 0 , (2.33) 

where 

14 = (cos0)I4+(sin 0)14. (2.34) 

As illustrated in Fig. p.l| , 6a is the component of the two-held perturbation vector 
along the direction of the background helds’ evolution. Conversely, huctuations 
orthogonal to the background classical trajectory represent non-adiabatic pertur¬ 
bations, and we dehne the “entropy held”, s, such that 

hs = (cos6*)5y — (sin6*)50. (2.35) 
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Figure 2.1: An illustration of the decomposition of an arbitrary perturbation into 
an adiabatic (^a) and entropy {6s) component. The angle of the tangent to the 
background trajectory is denoted by 9. The usual perturbation decomposition, 
along the 0 and y axes, is also shown. 


From this dehnition, it follows that s =constant along the classical trajectory, 
and hence entropy perturbations are automatically gauge-invariant |^. Pertur¬ 
bations in 6a, with 6s = 0, describe adiabatic held perturbations, and this is why 
we refer to a as the “adiabatic held”. 

The total momentum of the two-held system, given by Eq. (|2.9|), is then 


S(l,i = 


xSx,i = 




(2.36) 


and the comoving curvature perturbation in Eq. (p.l7|) is given by 




ij + H 


(j)6(j) + x6x 

02 + ;,^2 


tp H-OCT . 

a 


(2.37) 


This expression, written in terms of the adiabatic held, a, is identical to that 
given in Eq. p.27] ) for a single held. 

We can also write Eq. (|2.37|) as 


71 = (cos^ 0)71^ + (sin^ 9)71^ , 


(2.38) 
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where we define the comoving curvature perturbation for each of the original 
fields as 


H H 

TZj = iIj + —Sifi = —Qi. 

Vi Vi 


(2.39) 


However, even fields with no explicit interaction will in general have non-zero 
intrinsic entropy perturbations on large scales in a multi-field system due to their 
gravitational interaction, so that TZi for each field is not conserved. Although the 
intrinsic entropy perturbation for each field is still of the form given by Eq. ( |2.28| ), 
it is no longer constrained by Eq. (|2.14|) to vanish as /c —> 0. This is in contrast to 
the case of non-interacting perfect fluids, where it is possible to define a constant 
curvature perturbation for each fluid on large scales m 


The comoving matter perturbation in Eq. ( p.lOj) can be written as 


Cm = d ( 5(7 — a A ) — aba + 2Vs5s , 


(2.40) 


which acquires an additional term, compared with the single-field case, due to 
the dependence of the potential upon s, where 


Vs = (cos9)V^ — (sin6*)V^. 


(2.41) 


The perturbed kinetic energy of s has no contribution to first-order as in the 
background solution s = 0, by definition. 

The total entropy perturbation, Eq. ( p.23|) , for the two fields can be written 
as 


5 = 


3cT2(3hf(T + 214 


X 


X K 


a 16a — a A ) — a6a 


+ 3Ha^e6s 


Combining Eqs. (|2.14|) , (|2.4(]|) and (|2.42|) , we can write 

fC 2Vs 


5=_^_4$ 

67rGcT^[3hfcT-|-214] 


3d2 


6s. 


(2.42) 


(2.43) 


Comparing this with the single-field result given in Eq. ( p.29| ), we see that the 
entropy perturbation on large scales is due solely to the relative entropy pertur¬ 
bation between the two fields, described by the entropy field 6s. 


The change in the comoving curvature perturbation is given by |^, ^ 

-H 


, 

Ha^ 2 


f6(p 6x\ d / 0^ 


X J dt \ 02 ^2 


(2.44) 
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which can be expressed neatly in terms of the new variables: 


■ Hk\ 

7^ = + —06s, 

Ha^ 


a 


where 


a 


( 2 . 45 ) 


( 2 . 46 ) 


The new source term on the right-hand-side of this equation, compared with the 
single-held case, Eq. (|2.3CI| ), is proportional to the relative entropy perturbation 
between the two helds, 6s. Clearly, there can be signihcant changes to TZ on 
large scales if the entropy perturbation is not suppressed and if the background 
solution follows a curved trajectory, i.e., 0 7^ 0, in held space [Q. This can then 
produce a change in the comoving curvature on arbitrarily large scales (i.e., even 
in the limit fc —>■ 0) 


2g 


Equations of motion for the adiabatic and entropy held perturbations can be 
derived from the perturbed scalar held equations (|2.5|), to give 

60' 3H6& i~2 1 


= -2V„A + a 


. p . 

A + 3'iP + —{PE -aB) 


+ 2{e6sy -2^96s, 

a 


(2.47) 


and 


6s 3H6s I — -|- Vss — 0 I (js 


= - 2 - 
& 


&(6a — &A) — a6a 


(2.48) 


where 


14^ = (sin^ 0)V^x + (sin ‘^0W4>x + > (2-49) 

14s = (sin^ 9)V^4, - (sin 20)14^ -f (cos^ 0)V^^ . (2.50) 

When 9 = 0, the adiabatic and entropy perturbations decouple. If we employ 
the slow-roll approximation for the background helds, (p ~ —V^/3H and x — 
—14/377, we obtain 0^0. This rehects the fact that the rate of change of 0 is 
slow - instantaneously it moves in an approximately straight line in held space. 
But the integrated change in 0 cannot in general be neglected. Even working 
within the slow-roll approximation, helds do not in general follow a straight line 
trajectory in held space. The equation of motion for 6a then reduces to that for 
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a single scalar field in a perturbed FRW spacetime, as given in Eq. ( p.25D , while 
the equation for 6 s is that for a scalar held perturbation in an unperturbed FRW 
spacetime. 

The only source term on the right-hand-side in Eq. (p.48|) for the entropy 
perturbation comes from the intrinsic entropy perturbation in the a-held (see 
Eq. p.28|) . From Eqs. ( p.l4|) and (|2.4CI|) we have 


&(6a — a A) — a6a = 2&96s 






(2.51) 


47rGa^ 

and hence we can rewrite the evolution equation for the entropy perturbation as 


6s 3H6s -|- ( —IT -f- + 30 ) 5s — -: 


9 


a 


27rGa^ 


T. 


(2.52) 


Note that this evolution equation is automatically gauge-invariant and holds in 
any gauge. On large scales the inhomogeneous source term becomes negligible, 
and we have a homogeneous second-order equation of motion for the entropy 
perturbation, decoupled from the adiabatic held and metric perturbations. If the 
initial entropy perturbation is zero on large scales, it will remain so. 

By contrast, we cannot neglect the metric back-reaction for the adiabatic held 
huctuations, or the source terms due to the entropy perturbations. Working in 
the spatially hat gauge, dehning 


a 


Qa = SaQ = 6 a+ —- 0 , 

li 


(2.53) 


and using 


a 


Aq = iirG—Q ,, 


(2.54) 


we can rewrite the equation of motion for the adiabatic held perturbation as 


Qcr + “iHQa + 


+ 14 . - 0 ' 


(c?a^\ 


H 


Qo 


= 2{96sy -2 \Yy + ^\96s. 

(j H 


(2.55) 


When 0 = 0, this reduces to the single-held equation of motion, but for a curved 
trajectory in held space, the entropy perturbation acts as an additional source 
term in the equation of motion for the adiabatic held perturbation, even on large 
scales. 

In order for small-scale quantum huctuations to produce large-scale (super- 
Hubble) perturbations during inhation, a held must be “light” (i.e., overdamped). 
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The effective mass for the entropy held in Eq. (p.52|) is /i^ = 14^ + 36^. For 


/if > the huctuations remain in the vacuum state and huctuations on large 
scales are strongly suppressed. The existence of large-scale entropy perturbations 
therefore requires 


/if = + 302 < -H^ 


(2.56) 


2.3 Application to entropy/adiabatic correlations 
from inflation 


Equations (|2.52|) and (p.55|) are the key equations which govern the evolution of 
the adiabatic and entropy perturbations in a two held system. Together with 
constraint equations (p.51 ) and ( 2.54 ) for the metric perturbations, they form 
a closed set of equations. They allow one to follow the ehect on the adiabatic 
curvature perturbation due to the presence of entropy perturbations, absent in 
the single held model. This in turn will allow us to study the resulting correlations 
between the spectra of adiabatic and entropy perturbations produced on large- 
scales due to quantum huctuations of the helds on small-scales during inhation. 

A useful approximation commonly made when studying held perturbations 
during inhation, is to split the evolution of a given mode into a sub-Hubble regime 
{k > aH), in which the Hubble expansion is neglected, and a super-Hubble regime 
{k < aH), in which gradient terms are dropped. 

If we assume that both helds 0 and y are light (i.e., overdamped) during 
inhation, then we can take the held huctuations to be in their Minkowski vacuum 
state on sub-Hubble scales. This gives their amplitudes at Hubble crossing {k = 
aH) as 


I k=aH 






ei{k) , 


(2.57) 


where I = 0, y, Hj. is the Hubble parameter when the mode crosses the Hubble 
radius (i.e., Hk = k/a), and and are independent Gaussian random variables 
satisfying 


{ei{k)) = 0 , {ei{k)e*j{k')) = 6ij6{k - k '), 


(2.58) 


with the angled brackets denoting ensemble averages. It follows from our dehni- 
tions of the entropy and adiabatic perturbations in Eqs. ( p.31| ) and (|2.35| ) that 
their distributions at Hubble crossing have the same form: 


^^\k=aH 


H. 




e„{k ), (5s| 


Hy 


k=aH 




es{k) , 


(2.59) 
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where and are Gaussian random variables obeying the same relations given 
in Eq. (p.58|) , with I, J = a, s. 

Super-Hubble modes are assumed to obey the equations of motion given in 
Eqs. (p.55|) and (|2.52|) , which we will write schematically as 


d%6s) = 0 , 


(2.60) 

(2.61) 


where 0”'(Qo-) and 0^(5s) are obtained by setting fc = 0 on the left-hand side 
of Eqs. (p.55|) and (|2.52|) respectively, and S'^{6s) is given by the right-hand side 
of Eq. (|2.55|) . As remarked before, there is no source term for 6 s appearing on 
the right-hand side of Eq. ( p.52|) once we neglect gradient terms. The general 
super-Hubble solution can thus be written as 


= A+U{t) + A_f_{t) + P{t), (2.62) 

(5s = B+g+{t) + B_g_{t), (2.63) 


where the real functions f± and g± are the growing/decaying modes of the homo¬ 
geneous equations, 0”'(/±) = 0 and 0^{g±) = 0, and P{t) is a particular integral 
of the full inhomogeneous equation ( p.dOj) . Note that the slow-roll growing-mode 
solution /+ oc a/H. 

Henceforth we shall consider only slow-roll inflation where the evolution can 
be approximated by hrst-order equations [dropping 5s and Q„ in Eqs. ( |2.52|) 
and (|2.55|) ], so that we have0 

~ Af{t) + P{t), (2.64) 

5s ~ Bg{t). (2.65) 


We can, without loss of generality, take f = 1 = g and P = 0 when k = aH, 
so that the amplitudes of the growing modes at Hubble-crossing are given by 
Eqs. (p.59|) as 


A(k) 


B(k) 


Hk 


\/2k^ 


es{k) 


( 2 . 66 ) 


From Eq. ( p.6(]| ), we see that the amplitude of the particular integral P{t) at later 
times will be correlated with the amplitude of the entropy perturbation, P, and 
we can write P{t) = BP{t), where P{t) is a real function independent of the 
random variables Co-, e^. 

In order to quantify the correlation, we dehne 


27r2 

{x{k)y*{k')) = — Cxy 5{k - k'). (2.67) 

^We note that in non-slow-roll scenarios the decaying modes may not be negligible on super- 
Hubble scales, which could affect the correlations between adiabatic and entropy perturbations. 
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The adiabatic and entropy power spectra are given by 




QaQa — 


'Pss = Cs 


ih 

27r 

ih 

27r 


|f| + |P^ 


l/l, 


while the dimensionless cross-correlation is given by 

Cq„5s gP 


v/^^|/2| + |P2| 


( 2 . 68 ) 

(2.69) 


(2.70) 


Note that the adiabatic power spectrnm at late times is always enhanced if it is 
conpled to entropy pertnrbations [i.e., P{t) ^ 0, in Eq. (|2.64|) ], as the entropy 
held hnctnations at Hnbble-crossing provide an nncorrelated extra sonrce. 

As an illnstration, we consider the correlations in the adiabatic and entropy 
pertnrbations at the start of the radiation era, prodnced after donble inhation, 


as stndied in Ref. 
massive scalar helds is: 


The donble-inhation potential for two non-interacting bnt 


V = 



( 2 . 71 ) 


Following [Q, it is possible to parametrise the backgronnd scalar held trajectory 


in polar coordinates when both helds are slow-rolling: 


X - 






(2.72) 


where N = — ln(a/aend) is the nnmber of e-folds nntil the end of inhation. The 
backgronnd trajectory can then be expressed as: 


JV ~ JVn —---- 

“ °(cosa)2^^/(«^-b ’ 


(2.73) 


where R = The scalar held position angle, a, can be related to the 

scalar held velocity angle, 9, which we used to dehne the adiabatic and entropy 
perturbations: 


tan 6* ~ 


3Ha 




(2.74) 


The scalar held y is assumed to decay into cold dark matter while the scalar 
held 0 decays into radiation. The entropy/isocurvature at the start of the 
radiation-dominated era is described by 


y _ ^ Pc 3 ^P'y 

'rad = T • 

Pc 4 


( 2 . 75 ) 
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In Ref. [Q, it is shown how the snper-Hubble pertnrbations in the radiation 
era can be determined in terms of the pertnrbations during the inflationary era. 
The fluctuations in both 0 and y helds can contribute to both the adiabatic 
and entropy perturbations. The adiabatic component comes directly from the 
comoving curvature perturbation, TZ, at the end of inflation, and is given by 


7^I.ad ^ -V47rG 







Hk [{sin ak)e^{k) + {cosak)e^{k)] 


(2.76) 


The isocurvature perturbation at the start of the radiation-dominated era is re¬ 
lated to the entropy perturbation between the two helds at the end of inhation 


c ^ 2 2 1 ¥ 

*S^rad — Q rr ( • 1 

^ \ x (b 


which yields 


(2.77) 


Srad — -V4:7rG\ -rjHk [R^secak + cosec es{k ), 


(2.78) 


and 


n 


rad 


X 


— Nk ^ tan ak sin ak 
~ V47r(T\/ -rrHk , X 

V k yT? 2 ^iai?afc + l 


R^ tan^ ak 


+ 1 


ea{k) 


1-R2 

R‘^ tan ak 


es{k) 


(2.79) 


The entropy perturbation during the radiation era only depends on the entropy 
perturbation at Hubble-crossing during the inhationary era, while the adiabatic 
perturbation during the radiation era depends on both the adiabatic and entropy 
perturbations at Hubble-crossing. This is consistent with equations (|2.52|) and 
(|2.47|) , showing that the entropy perturbation sources the adiabatic perturbation 
on super-Hubble scales, but not vice versa. 

As both equations (|2.78|) and (p.79|) depend on the random variable e^, the 
adiabatic and entropy perturbations will be correlated, and we hnd 


C 


^rad 'S'rad 


(i?^ — 1 ) sin 2 ak 


2 SiU^ ttfc + COS^ Ofe) 


(2.80) 


This correlation is investigated fully in in terms of the usual scalar held per¬ 
turbation variables. An interesting point that can easily be seen from Eq. (|2.79|) 
is that 77rad will depend only on e„ ii R = m^/m^ = 1. Thus, there will be no 
correlation if 7? = 1. As can be seen from Eq. (p.73 ), a will be constant for R = 1 
and thus so will 0 ; a straight-line background trajectory will be obtained for 


R= 1. This is consistent with Eq. (|2.47|) , where it can be seen that the entropy 
component only sources the adiabatic component on large scales if 6 * 7 ^ 0 . 
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2.4 Conclusions 


We have introduced a new formalism in which to follow the evolution of adiabatic 
and entropy perturbations during inflation with multiple scalar helds. We decom¬ 
pose arbitrary held perturbations into a component parallel to the background 
solution in held space, termed the adiabatic perturbation, and a component or¬ 
thogonal to the trajectory, termed the entropy perturbation. We have rederived 
the held equations in terms of these rotated helds in Eqs. ( p.52|) and (|2.55|) . These 
show that the adiabatic perturbation on large scales can be driven by the entropy 
perturbation, while the entropy perturbation itself obeys a homogeneous second- 
order equation on super-Hubble scales. There can only be signihcant change 
in the large-scale comoving curvature perturbation if there is a non-negligible 
entropy perturbation, and if the background trajectory in held space is curved. 

Our formalism can be applied to evaluate the correlation between the adia¬ 
batic and entropy perturbations at the end of inhation. As an example we con¬ 
sidered the example of two non-interacting helds in double inhation, calculating 
the cross-correlation between the adiabatic and entropy perturbations. 
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Chapter 3 
Preheating 


3.1 Introduction 


Standard inflationary models must end with a phase of reheating during which 
the inflaton, 0, transfers its energy to other fields. Reheating itself may begin 
with a violently nonequilibrium “preheating” era, when coherent inflaton oscil¬ 
lations lead to resonant particle production (see and refs, therein). Until 


recently, preheating studies implicitly assumed that preheating proceeds with¬ 
out affecting the spacetime metric. In particular, causality was thought to be a 
“silver bullet,” ensuring that on cosmologically relevant scales, the non-adiabatic 
effects of preheating could be ignored. 

During preheating metric perturbations may be resonantly amplified on all 
length scales Causality is not violated precisely because of the 

huge coherence scale of the inflaton immediately after inflation ^ (see also 
H). Strong preheating (with resonance parameter g 3> 1; see Chapter |l] for 
overviews and notation) typically leads to resonant amplification of scalar metric 
perturbation modes <hfc, possibly even on super-Hubble scales (i.e., k/aH 1). 

Understanding when super-Hubble scales are amplified is crucial since pre¬ 
heating can lead to e anisotropies in the CMB. Observational limits rule out 
those models that produce unbridled nonlinear growth, but models which pass the 
metric preheating test on COBE scales may nevertheless leave a non-adiabatic 
signature of preheating in the CMB. 


3.2 Entropy suppression 

In this section we use the entropy/adiabatic decomposition of the perturbation 
equations developed in Chapter |] to investigate the dynamics of super-Hubble 
perturbations during a period of preheating at the end of inflation. We consider 
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three models, encompassed by the general effective potential 


1/ = -m" 
2 


A 

+ 4^ 


1 


i z I ~2 i3 

X + 9 X ■ 


(3.1) 


The essence of preheating lies in the parametric amplification of field perturba¬ 
tions due to the time-dependence of their effective mass, e.g., 

In the simplest cases, the inflaton 0 simply oscillates at the end of inflation. 

Preheating typically amplihes long-wavelength modes preferentially. As dis¬ 
cussed in amplihcation of super-Hubble modes does not lead to a 

violation of causality, due to the super-Hubble coherence of the inflaton oscilla¬ 
tions set up by the prior inflationary phase. If TZ is amplified on super-Hubble 
scales, this will alter the resulting imprint on the anisotropies of the CMB, and 
break the simple link between CMB observations and inflationary models. 

We consider first the case where the inflaton is massive (m ^ 0) and neglect 
its self-interaction (A = 0). The traditional resonance parameter for the strength 
of preheating at the end of inflation is 


q. = 


9 ^ 

Am? ’ 


(3.2) 


where 0o is the initial value of 0 at the beginning of preheating. In the massive 
case, where modes move through the resonance bands of the Mathieu chart, 
and for inflation at high energies where the expansion of the universe is very 
vigorous, q needs to be much larger than one if the parametric resonance is to 
be efficient [|^. It is possible to have large q even for small coupling, -C 1, as 
m <C 00 ~ Mpi. We can write the effective mass of the y during inflation as 




3q M] 




TT 0^ 


(3.3) 


It then follows from Eq. (|3.3|) that y must be heavy during inflation for this 
simple potential if efficient preheating is to be obtained. 

Any change in the curvature perturbation IZ on very large scales must be due 
to the presence of non-adiabatic perturbations. In Q, it was shown how, 
if during inflation with A = 0 = ^, then the y held and hence any 

non-adiabatic perturbations on large scales are exponentially suppressed during 
inhation, and no change to IZ occurs before backreaction ends the resonance. 

However, when ^ 7 ^ 0, the y held will have a nonzero vacuum expectation 
value (vev) during inhation even along the valley of the potential. In the slow-roll 
limit for 0, this vev is determined by = 0, which gives (see Eq. 2.2 ) 


9 

X ~- 7 ' 

9^ 


(3.4) 
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The g coupling has the effect of rotating the valley of the potential - which the 
attractor trajectory approximately follows - from y = 0 , through an angle 


9^ 


(3.5) 


where, to ensure that the chaotic inflation scenario is not drastically altered, we 
assume 0 


9 


< 1 . 


(3.6) 


The effect of g is to change the attractor for both y and 5y during inflation, 
since the y and 6 x equations of motion (Eqs. p.2| and |2.5| ) gain inhomogeneous 
driving terms proportional to g‘^(j)^. This does not necessarily imply that TZ will 
be amplihed by preheating at the end of inflation as purely adiabatic perturba¬ 
tions along the slow-roll attractor now have a component along y as well as 0. In 
order to determine whether or not the evolution of the comoving curvature per¬ 
turbation, 71, on super-Hubble scales is affected, we need to follow the evolution 
of the entropy field perturbatiorQ, dehned by Eq. (p.35|) , which gives 


5s ~ (5y -1- ^S(j). 
9 ^ 


(3.7) 


In the limit 0 we recover 5s —>■ 5y. Crucially, the evolution equation ( |2.52| ) 

for the entropy perturbation has no inhomogeneous terms in the long-wavelength 
{k —> 0 ) limit, even for ^ 7 ^ 0 , and entropy perturbations will only be non- 
negligible on super-Hubble scales if the entropy held is light during inhation. 

In the slow-roll limit and on large scales, the evolution equation (|2.52|) for the 


entropy perturbation has the approximate solution [39 


where 


Ss (X a i 

an 


.2 = 9 _ 

4 ’ 


(3.8) 


(3.9) 


and the effective mass of the entropy held, is dehned in Eq. (|2.56|) . The power 
spectrum of entropy perturbations is 




3-2Re(i/) 


(3.10) 


^ From Eq. ( ^.45 ) we see that 95s must be non-zero to change TZ on large scales. Because 
0 ~ 0, from Eq. (3.5), the entropy remains decoupled from the adiabatic perturbation during 
slow-roll inflation in this model. But at the end of inflation, during preheating, 0 7^ 0. 




















The real part of z/ vanishes for >9/4, leaving a steep bine spectrnm, 

which is exponentially snppressed with time. 

Using Eqs. CT , (O), (ED, and ( p.6|) , one hnds that 
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3qMpy 

TT(Pl 


(3.11) 


has a local minimnm for ^ = 0. Thns the additional g term in Eq. (|3.1|) 
serves to increase the entropy mass relative to the Hnbble parameter, and so does 
not avoid the snppression of the entropy pertnrbation. The g term therefore does 
not signihcantly alter the spectral index of the spectrnm of entropy pertnrbations, 
which remains steep if g >> 1. The strongly bine spectrnm implies that non-linear 
backreaction is dominated by small-scale modes, which go nonlinear long before 
the cosmological modes, implying that resonance ends before 71 changes []^ 

We have also integrated the held eqnations nnmerically to avoid relying on 
any slow-roll-type approximations. To nnmerically evalnate the entropy pertnr¬ 
bation, one conld simnlate the original pertnrbation variables 50 and 6 x, nsing 
Eq. (|2.5|) , and then work ont 6 s algebraically via Eq. ( p.35|) . However, this 
approach is prone to nnmerical instability when the entropy perturbation is sup¬ 
pressed. To illustrate this, we take g = 8 x and q = 3.8 x 10^ After about 

60 e-folds of inhation, one can see analytically that 6 s ~ 10 “^°. Numerically, 
6 xcos6 ~ 50 sin 6* ~ 10 “® during inhation. So in order to obtain a high enough 
accuracy to model the suppression of 5s, we require that 5% cos 6* and 50 sin 6* 
have to be simulated to a relative accuracy of 10 ®/10 = 10 This means 

approximately 32 signihcant hgures are needed, which is beyond the capability 
of standard numerical ordinary diherential equation integration routines. 

If instead we use the new adiabatic and entropy held perturbations and inte¬ 
grate Eqs. ( |2.52| ) and (|2.55|) , then this numerical instability does not occur, since 
one no longer needs to hnd the diherence between two nearly equal quantities. 
Simulation results using these equations are compared with the results using the 
old held perturbation equations 


in Fig. |3]^. The simulations show that the 
growth in 71 is driven by 5s, in concordance with Eq. ( ^.45| ). As can be seen, the 
numerical result using the held perturbation equations fails to track the expo¬ 
nential decay of the entropy during inhation and thus underestimates the delay 
in the growth of 7Z. 

In practice, we hnd a similar instability if we try to construct the gauge- 
invariant metric perturbation, T, required in Eq. ( |2.52|) in terms of the constraint 
Eq. (|2.51|) . This includes the intrinsic entropy perturbation in the a held, which 
does become small at late times/large scales, but results from the diminishing 
diherence between hnite terms. It is more stable numerically to follow the value 


29 




















Figure 3.1: Numerical simulations of the entropy and comoving curvature pertur¬ 
bations during inflation and preheating, with A = 0, (7 = 2 X 10-3, g = 8 x 
and m = 10-®Mpi. The ‘new’ prehx indicates that the field perturbations were 
evaluated by numerically integrating Eqs. (|2.52|) , and ( p.55|) , while the ‘old’ pre¬ 
fix indicates that the perturbations were evaluated by integrating the original 


field equations (p.5|) . We have not included any higher-order corrections such as 
backreaction from small-scale perturbations which would shut down the resonant 
amplification of 6 s at some point. 


of T at late times using the evolution equation 

+ I j = inGaQ ^, 


(3.12) 


which can be obtained from the definition of T given in Eq. (p.l4|) and the metric 
constraint equations (|2.7|) and (p.ll|) . 

Note that the adiabatic/entropy decomposition becomes ill-defined if d = 0, 
i.e. both helds stop rolling, and this can cause numerical instability during pre¬ 
heating if the trajectory is conhned to a narrow valley. This can occur, for in¬ 
stance, when ^ = 0 and only the (j) field oscillates. The original held perturbations 
6 (j) and 6 x remain well-dehned, although the comoving curvature perturbation TZ, 
dehned in Eq. (|2.37 ) becomes singular when a = 0 [^. This does not happen 
for the simulation results shown in Fig. with g ^ 0 where the helds oscillate 
in a two-dimensional potential well. 

The massive inhaton potential (m ^ 0) safeguards the conservation of 77. by a 
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bootstrap effect: if preheating is strong, g S> 1, then the entropy perturbation is 
heavy during inflation; on the other hand, if the entropy is light during inflation, 
then g < 1 and preheating is very weak. This is not altered by a rotation of 
the traiectory in held space (g ^ 0) as can most quickly be seen by noting, from 
Eqs. and (|^), that 

Kro- + 14s = V(i><t> + Vxx ■ (3.13) 

Thus if the y held is very massive S> ff^), we must have 14-o- + 14s S> 

For slow-roll inhation we require 14^cr and hence 14s 3> H^. 


3.3 Unsuppressed entropy 


The suppression does not necessarily occur in massless (m = 0) self-interacting 
(A 7 ^ 0) inhation models |^, This latter class of models is almost con¬ 
formally invariant, allowing analytical results from Floquet theory to be applied. 
The Floquet index, which determines the rate of exponential growth, can 
reach its maximum as k/aH 0, when g^/X = 2n^ for integer n, thereby im¬ 
plying maximum growth for the longest-wavelength perturbations. Assuming 
slow-roll inhation driven by F ~ A0^/4, we see from Fq. (|3.13|) that 14o- -l- 14s > 
14 x = and thus that the entropy held is massive ( 14 s > 9if^/4) whenever 


9 ^ 

^ > Svr 
A 


MJ, 


(3.14) 


However, we can have resonance at large scales for n = 1 and g^ j A = 2, when the 
entropy held need not be heavy during inhation and no exponential suppression 
takes place, so that the subsequent growth of 7^ is explosive [Q. The growth 
of 'R. occurs before backreaction can shut oh the resonant growth of the entropy 
perturbations bs 0, §g, 0,0. Although the region of parameter space around 
g^jX = 2 is thus ruled out, the same does not hold for g^jX S> 1, since the 
entropy held is then heavy during inhation and 5s is again suppressed. 

Models where the entropy ehective mass is simply very small during inhation 
but then becomes large at preheating can also ehect large scales. 


3.4 New cosmological effects 


metric preheating can lead to a host of 


Beyond the ehects discussed in [ 
interesting new ehects. 

• The growth of implies amplihcation of isocurvature modes in unison with 
adiabatic scalar modes on super-Hubble scales. Preheating thus yields the pos¬ 
sibility of inducing a post-inhat ionary universe with both isocurvature and adi¬ 
abatic modes on large scales. The ehect of having mixtures of isocurvature and 
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adiabatic perturbations on the CMB is discussed in Chapter 
• Because the metric perturbations can go nonlinear, whether on sub- or super- 
Hubble scales, the corresponding x density perturbations 6 typically have non- 
Gaussian statistics. This is simply a reflection of the fact that — 1 < 5 < cxd, 
so that the distribution of necessity becomes skewed and non-Gaussian. Fur¬ 
ther, when (x) = 0 during inflation, x perturbations in the energy density will 
necessarily be non-Gaussian (chi-squared distributed), even if 6 xk is Gaussian 
distributed, since stress-energy components are quadratic in the fluctuations (see 

e.g. IZl). 


• Another new feature we can identify is the breaking of conformal invariance. 
Once metric perturbations become large on some scale, the metric on that scale 
cannot be thought of as taking the simple Friedmann-Robertson-Walker (FRW) 
form, and conformal invariance is lost. This is particularly important for the 
production of primordial magnetic fields, which are usually strongly suppressed 
due to the conformal invariance of the Maxwell equations in a FRW background. 
The coherent oscillations of the inflaton during preheating further provide a nat¬ 
ural cradle for producing a primordial seed for the observed large-scale magnetic 
fields. A charged inflaton field, with kinetic term , will couple to 

electromagnetism through the gauge covariant derivative This 

will naturally lead to parametric resonant amplification of the existing magnetic 
held, which could produce large-scale coherent seed fields on the required super- 
Hubble scales without fine-tuning |^. (Note that a tiny seed held must exist 


during inflation due to the conformal trace anomaly and one-loop QED correc¬ 
tions in curved spacetime [ff^].) 


3.5 Conclusion 

The effect of preheating on the large-scale curvature perturbation has been ex¬ 
amined using the formalism developed in Ghapter The mass of the entropy 
held during inflation is a crucial quantity. If the entropy held is heavy, then any 
fluctuations on large scales are suppressed to negligible values at the beginning 
of preheating. This squeezing of the entropy perturbation is most accurately 
modelled numerically using our evolution equation for the entropy perturbation. 
If it is estimated from the usual held equations, it may contain large numerical 
errors when there is a non-trivial background trajectory in held space. 

For models where efficient preheating can occur with a light entropy dur¬ 
ing inflation, large scale perturbations are effected by preheating. A range of 
cosmological implications for such models was discussed. 
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Chapter 4 

Correlated adiabatic and entropy 
perturbations and the Cosmic 
Microwave Background 


4.1 Introduction 


Increasingly accnrate measnrements of temperatnre anisotropies in the cosmic 
microwave background sky offer the prospect of precise determinations of both 
cosmological parameters and the nature of the primordial perturbation spectra. 
The recent Boomerang m , DASI [0 and Maxima IJT^ data have shown evidence 
for three peaks in the CMB temperature anisotropy power spectrum as expected 
in inflationary scenarios. In this context the CMB data support the current ‘con¬ 
cordance’ model based on a spatially flat Friedmann-Robertson-Walker universe 
dominated by cold dark matter and a cosmological constant [^. In addition, 
the CMB data no longer shows any signs of being in conflict with the big bang 


nucleosynthesis data |76 


In the studies which have estimated the cosmological and primordial param¬ 
eters with these new data sets, only the case of purely adiabatic perturbations 
has been considered so far, i.e. the perturbation in the relative number densities, 
6 n/n, of different particle species is taken to be zero. Although this assumption 
is justihed for perturbations originating from single held inhationary models, it 
does not necessarily follow when there is more than one held present during in- 
hation (see for example ^ |^). Other possible primordial modes are 

isocurvature P^ , P6[| (also referred to as “entropy”) modes in which the particle 
ratios are perturbed but the total energy density is unperturbed in the comoving 
gauge. 

Most previous studies have examined the extent to which a statistically inde¬ 
pendent isocurvature contribution to the primordial perturbations may be con¬ 
strained by CMB and large-scale structure data [03, It has recently been 


33 


















shown that multi-held inhationary models in general produce correlated adia¬ 
batic and isocurvature perturbations m 0 ^ Hi- These correlations can dra¬ 
matically change the observational ehect of adding isocurvature perturbations 


. Up until now, only the case of scale-invariant correlated adiabatic 
and entropy perturbations has been considered. Trotta et al. found (with an 
earlier CMB dataset) that in this case the cold dark matter (CDM) isocurvature 
mode was likely to be very small if not entirely absent, though they did hnd 
that a neutrino isocurvature mode contribution M, M was not ruled out. In 


this Chapter we examine whether a correlated CDM isocurvature mode is better 
favoured by the recent CMB data when a tilted power law spectrum is allowed. 


4.2 Theory 


Non-adiabatic perturbations are produced during a period of slow-roll inhation 
in the presence of two or more light scalar helds, whose ehective masses are less 
than the Hubble rate. On sub-horizon scales, fluctuations remain in their vacuum 
state so that when fluctuations reach the horizon scale their amplitude is given 
by ~ (hf*/27r)ai where the subscript * denotes horizon-crossing and a* are 
independent normalised Gaussian random variables, obeying { 0 ,^ 0 ,j) = 6ij. The 
total comoving curvature and entropy perturbation at any time during two-held 
inhation can quite generally be given in terms of the held perturbations, along 
and orthogonal to the background trajectory, (see Chapter H) 

iZ oc cos 6* 501 -|-sin 6* 502 ) (4-1) 

S oc — sin 0 501-t-cos 0 502 , (4.2) 


where 9 is the angle of the inhaton trajectory in held space. Although the curva¬ 
ture and entropy perturbations are uncorrelated at horizon-crossing, any change 
in the angle of the trajectory, 9, will begin to introduce correlations. Further 
correlations may be introduced by the model dependent dynamics when inhation 
ends and the helds’ energy is transformed into radiation and/or dark matter. 
The comoving curvature perturbation, T^rad, on large-scales during the radiation- 
dominated era is related to the conformal Newtonian metric perturbation, $, by 
’^rad = 3<F/2. The isocurvature perturbation is 

(o _ 5pcdm 3 5 p.y 

‘^rad ~7 

Pcdm 4 


and remains constant on large scales until it re-enters the horizon. On large scales 
the CMB temperature perturbation can be expressed in terms of the primordial 
perturbations 


5T 

~T 


1 

5 



(4.3) 
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The general transformation of linear curvature and entropy perturbations 
from horizon-crossing during inflation to the beginning of the radiation era will 
be of the form 


^rad 

‘^rad 


1 Ttis 
0 Tss 


5* 


(4.4) 


Two of the matrix coefficients, = 1 and = 0, are determined by the phys¬ 
ical requirement that the curvature perturbation is conserved for purely adiabatic 
perturbations and that adiabatic perturbations cannot source entropy perturba¬ 
tions on large scales [^. The remaining terms will be model dependent. If the 
fields and their decay products completely thermalize after inflation then Tss = 0 
and there can be no entropy perturbation if all species are in thermal equilibrium 
characterised by a single temperature, T. This means that it is unlikely that 
a neutrino isocurvature perturbation could be produced by inflation unless the 
reheat temperature is close to that at neutrino decoupling shortly before primor¬ 
dial nucleosynthesis takes place. On the other hand, a cold dark matter species 
could remain decoupled at temperatures close to, or above, the supersymmetry 
breaking scale yielding Tss- The simplest assumption being that one of the fields 
can itself be identified with the cold dark matter ^ . 

The slow evolution (relative to the Hubble rate) of light fields after horizon¬ 
crossing translates into a weak scale dependence of both the initial amplitude of 
the perturbations at horizon crossing, and the transfer coefficients T-jis and Tss- 
Parameterising each of these by simple power-laws over the scales of interest, 
requires three power-laws to describe the scale-dependence in the most general 
adiabatic and isocurvature perturbations. 


"^rad — T -A-gk^^CLs 1 

5rad = . 


(4.5) 

(4.6) 


The generic power-law spectrum of adiabatic perturbations from single field in¬ 
flation can be described by two parameters, the amplitude and tilt, A and n. Un¬ 
correlated isocurvature perturbations require a further two parameters, whereas 
we now have in general six parameters. The dimensionless cross-correlation 

^ _ (^rad^rad) _ sign(H) Agk^^ 
is in general scale-dependent. 


4.3 Likelihood analysis 

We will investigate in this Chapter the restricted case where all the spectra share 
the same spectral index and hence A is scale-independent. This might naturally 
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arise in the case of almost massless fields where the scale-dependence of the 
held perturbations is primarily due to the decrease of the Hubble rate during 
inhation, which is common to both perturbations and yields n* < 0. In the 
following analysis we also allow n* > 0, but we shall see that blue power spectra 
of this type are not favoured by the data. 

We then have four parameters, A = B, A, and n describing the 

effect of correlated perturbations, where n = 1 -|- 2 ni is dehned to coincide with 
the standard dehnition of the spectral index for adiabatic perturbations. We 
leave an investigation of the full six parameters for future work. 

By dehning the entropy-to-adiabatic ratio B* = B/A the parameter A be¬ 
comes an overall amplitude that can be marginalized analytically (see below). In 
the following, to simplify notation, we write A = 1 and drop the star from B*. 
We limit the analysis to H > 0 and 0 < A < tt, since there is complete symmetry 
under A ——A and under {B —B,A —>■ vr — A). Further, we allow three 
background cosmological parameters to vary, lji, = , uic = Vtcdmh’^ , and 

Ha where ^lb,cdm,A is the density parameter for baryons, CDM and the cosmolog¬ 
ical constant, respectively. Since we assume spatial flatness, the Hubble constant 
is 

2 a;c -f Wb 

Our aim is therefore to constrain the six parameters 

= {H, A, 77., (Ub, (Uc, Haj" , 


by comparison with CMB observations. We consider the COBE data analysed 
, and the recent high-resolution Boomerang [Q and Maxima data [|T^ 


m 


In 

order to concentrate on the role of the primordial spectra (and limit the numer¬ 
ical computation required) we will fix the reionisation history (no reionisation), 
neutrino masses (zero) and spatial curvature (zero). We will also neglect any 
contribution from tensor (gravitational wave) perturbations. 

We will use a CMBFAST code modified in order to allow correlated 
perturbations to calculate the expected CMB angular power spectrum, Ci, for 
all parameter values. (Our Ci is dehned as Ci = l{l + l)Cf/(27r) where 
is the square of the multipole amplitude). The computations required can be 
considerably reduced by expressing the spectrum for a generic value of B and A 
as a function of the spectra for other values. Let ns denote the purely adiabatic 
and isocurvature spectra when H = 1 as [Qjad and [Qliso respectively, and the 
correlation term for totally correlated perturbations H = 1,A = 0 as [C/jcorr • 
Then we can write the generic spectrum for arbitrary B and A as 


Q = IQU + B"|G]i.o + 2Bcos A[C,|< 


( 4 . 8 ) 


We can obtain [C/jcorr from Eq. (|4.8|) and using any HcosA ^ 0. The library 
spectra [C/jad and [C/jiso and [C^jcorr can then be used to evaluate Ci for any B and 
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A. A different set of library spectra will be needed for each set of cosmological 
parameters. When rii 7 ^ then A is not generally scale independent and so it 
would be necessary to evaluate the shape of the cross-correlation spectra [Gjcorr 
for each form of A(fc), but one can always perform the scaling with respect to B 
analytically. 

The remaining input parameters requested by the CMBFAST code are set 
as follows: Tcmb = 2.72QK, Y^e = 0.24, = 3.04, Tc = 0. In the analysis 

of 0 A, the optical depth to Thomson scattering, was also included in the 
general likelihood and, in the flat case, was found to be compatible with zero at 
slightly more than Icr. Therefore here, to further reduce the parameter space, 
we assume Tc vanishes. We did not include the cross-correlation between band 
powers because it is not available, but it should be less than 10 % according to 
An offset log-normal approximation to the band-power likelihood has been 


12 


advocated by and adopted by |^, |^, but the quantities necessary for its 
evaluation are not available. Since the offset log-normal reduces to a log-normal 
in the limit of small noise we evaluated the log-normal likelihood 


-2\og L{aj) = 






cr; 


( 4 . 9 ) 


where Ze = logC^, the subscripts t and d refer to the theoretical quantity and to 
the real data, Ce are the spectra binned over some interval of multipoles centered 
on £i, Gi are the experimental errors on and the parameters are denoted 
collectively as 

The overall amplitude parameter A can be integrated out analytically using 
a logarithmic measure d\ogA in the likelihood. Analogously, an analytic inte¬ 
gration can get rid of the calibration uncertainty of the Boomerang and Maxima 
data (see BH). to obtain the hnal likelihood function that we discuss in the 
following. We neglected beam and pointing errors, but we checked that the re¬ 
sults do not change signihcantly even increasing the calibration errors by 50%. 
We assume a linear integration measure for all the other parameters. 

In order to compare with the Boomerang and Maxima analyses we assume 
uniform priors as in W1 


with the parameters conhned in the range 


Bg( 0,3), AG(0,7r), ne (0.6,1.4) 

cUb G (0.0025, 0.08), G (0.05, 0.4), IIa e (0, 0.9). 


As extra priors, the value of h is conhned in the range (0.45, 0.9) and the universe 
age is limited to > 10 Gyr as in |]^. A grid of ~ 10, 000 multipole CMB spectra is 


used as a database over which we interpolate to produce the likelihood function. 

Figure 1 shows the one of the best cases in our database, corresponding to 
(B, n, (Ufe, cUc, Ga) = (0.3,0.8,0.02,0.1,0.7) and A = 0. The case A = 7 r /4 pro¬ 
vided an equally good ht. In the hgure the adiabatic ([Qjad), entropy {B^[Ci\iso) 
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Figure 4.1: A decomposition of the best fit spectra (solid line) which has A = 0. 
The adiabatic (dashes), entropy (dash-dot) and correlation (dots) contributions 
are also plotted. 


and correlated (25 cos A[Q]corr) components are shown. As can be seen, the ef¬ 
fect of adding a positively correlated component is to reduce the height of the 
low-/ plateau relative to the acoustic peaks . This is in contrast to the uncor¬ 
related case where the addition of entropy perturbations reduces the peak height 
relative to the plateau. 

We found a near-degeneracy between B and n when A = 0: the effect of 
adding maximally correlated isocurvature perturbations mimics an increase in the 
primordial slope. This makes clear the importance of varying n when studying 
correlated isocurvature perturbations: a lower n allows a larger B to be consistent 
with the CMB data. 

In Fig. 2 we plot a series of two-dimensional likelihood functions; all the other 
parameters have been marginalized in turn. The contour lines of the cosmological 
parameters ujb and ujc are almost parallel to B for 5 < 1. This means that the 
isocurvature perturbations do not alter signihcantly the best estimates for these 
cosmological parameters. On the other hand, increasing B moves the region of 
conhdence for Oa and of n toward smaller values. 

In Fig. 3 we plot the one-dimensional likelihood functions obtained by marginal¬ 
izing all the remaining parameters. Panel a shows that the contribution of isocur¬ 
vature perturbations can be as large as the adiabatic perturbations, or even larger: 
we hnd that B < 1.5 to 95% c.h. In contrast, if the isocurvature perturbations 
are uncorrelated, their fraction cannot exceed 70% {B < 0.7) to the same c.h. 
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Figure 4.2: Two-dimensional likelihood contour plots. The contours enclose 
40%, 86% and 99% of the likelihood and the stars mark the peaks. See text for 
explanation. 


It is intriguing to observe that the likelihood of B peaks around 0.3: that is, a 
non-zero contribution of isocurvature perturbations is more likely than a vanish¬ 
ing contribution. In the same panel we show as a long dashed line the likelihood 
assuming experimental errors reduced to one third, a precision within reach of 
the forthcoming satellite experiments: the curve shows that this level of precision 
would allow the detection of a hnite isocurvature contribution. Equally interest¬ 
ing, in panel b we see that the likelihood of the correlation cosA peaks near 
unity (maximal correlation), but has a not negligible probability everywhere in 
its domain. The likelihood functions for n and flA move toward smaller values, 
as anticipated, while the CDM and the baryon density estimates remain largely 
unaffected. The average values are 

n = 0.87 ±0.1, (Ub = 0.023 ± 0.005, 

= 0.12 ± 0.04, IIa = 0.63 ± 0.13. 
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Figure 4.3: One-dimensional likelihood functions in arbitrary units. Green 
(light) solid lines for the purely adiabatic models ( B = 0); blue short-dashed 
lines for uncorrelated fluctuations (cos A = 0); red (dark) solid lines for correlated 
fluctuations. See text for further explanation. 


4.4 Conclusions 


By contrast, Enqvist et al found that a large uncorrelated isocurvature con¬ 
tribution is only consistent with blue tilted slopes. The reason for this difference 
is that correlations can cause the acoustic peak height to increase relative to the 
Sachs Wolfe plateau (see Fig. 1) unlike the case of independent perturbations 
where the relative height always decreases. Trotta et al found that the CMB 
data was not consistent with a signihcant CDM isocurvature contribution because 
they restricted the primordial slope, n, to be unity. As can be seen from Fig. 2 
our n = 1 likelihood contours also indicate a very low isocurvature contribution. 
But when n is allowed to be less than one the isocurvature contribution can be 
even larger than the adiabatic contribution. 

As can be seen from Figs. 2 and 3 our estimates of ujb and ujc are virtually 
unaffected by the addition of correlated CDM isocurvature perturbations. Thus, 
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in our model, the nature of the isocurvature component can be investigated almost 
independently of the composition of the matter component. 

The main conclusion of this Chapter is that Boomerang and Maxima are 
consistent with a large correlated CDM isocurvature perturbation contribution 
when the spectral slopes are tilted to the red (n < 1). The higher precision of 
future satellite data has the potential to detect the isocurvature contribution, if 
any, thereby showing that inflation was not a single-held process. 
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Chapter 5 

Non-adiabatic perturbations 
from brane world effects 


5.1 Introduction 


According to string and M-theory, gravity is a higher-dimensional theory, reduc¬ 
ing to Einstein’s four-dimensional theory of general relativity at low enough ener¬ 
gies. In the brane-world scenario, the standard model matter fields are confined 
to a 3-brane in 1 -|- 3 -|- d dimensions, while the gravitational field can propagate 
in the bulk, i.e., also in the d extra dimensions, being localized at the brane 
at low energies. Recent developments show that the d extra space dimensions 
need not be small, or even compact, thus allowing the intriguing possibility that 
corrections could occur even at TeV scales. 

These exciting theoretical developments may offer a promising route towards 
a quantum gravity theory. However, as well as theoretical elegance, they must 
also pass the increasingly stringent tests provided by cosmological observations. 
Primarily, this involves developing higher-dimensional perturbation theory and 
then applying it to analyze the generation and evolution of density and tensor 
perturbations on the brane, leading to a prediction of the CMB anisotropies and 
galaxy distribution. 

This is an ambitious and difficult programme, but initial steps have already 
been taken, at least in the case of a particular class of models that generalize the 
Randall-Sundrum models [^. Large-scale adiabatic density perturbations from 
inflation on the brane have been computed (see also [^), using the con¬ 
servation of the curvature perturbation on uniform-density hypersurfaces. This 
conservation follows from adiabaticity and the conservation of energy-momentum 
on the brane, and is independent of the form of the field equations 


In 


the 


backreaction effect of metric fluctuations in the hfth dimension was neglected. In 
the general case, i.e., incorporating also the fluctuations in the nonlocal quantities 
that carry the bulk influence onto the brane, it has been shown that large-scale 
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density perturbations contain a closed system on the brane—and thus can in prin¬ 
ciple be evaluated purely from initial conditions on the brane, without knowledge 
8 ^]. Note that not all large-scale scalar perturbations can be 
the relation between the metric potentials A and 'h is 


of bulk dynamics 


computed intrinsically 
mediated by anisotropic stress imprinted on the brane by the bulk Weyl tensor 
and this stress cannot be evaluated without solving the bulk equations |^. In 
this Chapter, we solve the closed system described in to hnd the evolution of 
large-scale density perturbations on the brane. We show that extra-dimensional 
effects introduce a non-adiabatic mode on the brane. 

A general perturbation formalism has been developed |^, encompassing 
equations on the brane and in the bulk, and in principle able to describe all 
scales. However, the general equations are extremely complicated, in particular 
since the mode equations are partial differential equations. A hrst application of 
the equations has been made to large-scale tensor perturbations from inflation 
on the brane [^]. Unlike the scalar case, large-scale tensor perturbations cannot 
be evaluated without the bulk perturbation equations. 

We develop the outline argument presented hrst in |^, and analyze large- 
scale density perturbations and their evolution, from after Hubble-crossing in 
inhation through the radiation era. This provides part of the information needed 
for predicting the large-angle scalar anisotropies generated in CMB temperature, 
and seeing how the bulk effects modify general relativistic predictions. How¬ 
ever, the Sachs-Wolfe (SW) effect cannot be computed without knowledge of the 
Weyl anisotropic stress [Q. It is possible to estimate the SW effect by mak¬ 
ing assumptions about the Weyl anisotropic stress but a complete solution 
requires solving the full 5D perturbation problem. 

We show that in general, the perturbation $ (a covariant analog of the 
Bardeen metric perturbation) is no longer constant during high-energy inha¬ 


tion, but grows. However, is constant during the radiation era, as in general 
relativity, except at most in the early radiation era, if the energy density is still 
high relative to the brane tension. 


5.2 Brane dynamics 


We follow the 1-1-3 covariant approach and notation of [^| which is diherent 
from (but ultimately equivalent to) the metric-based approach to perturbations 
described in Chapter |^. The 5-dimensional (bulk) held equations are 


Gab — 


—^Qab + <5(x) {—Xqab + Tab} 


(5.1) 


where tildes denote th^bulk generalization of standard general relativity quan¬ 
tities, and 'k^ = Svr/Mp, where Mp is the fundamental 5-dimensional Planck 
mass, which is typically much less than the ehective Planck mass on the brane. 
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Mp = 1.2 X 10^® GeV. The brane is given by x = 0, so that a natural choice 
of coordinates is where x^ = are spacetime coordinates on 

the brane and the brane a moving hypersurface. The brane tension is A, and 
dAB = 'qab — nAUB is the induced metric on the brane, with ua the spacelike 
unit normal to the brane. Standard-model matter helds conhned to the brane 
make up the brane energy-momentum tensor Tab (with Tab^^ = 0). The bulk 
cosmological constant A is negative, and is the only 5-dimensional stress energy. 
(See [Q for the modihcation of this approach in the case where there is also a 
scalar held in the bulk.) 

The most general background bulk with homogeneous and isotropic induced 
metric on the brane is the 5-dimensional Schwarzschild-anti-de Sitter metric [PT[| , 
with the black hole mass leading to an effective radiation-like correction to the 
Friedmann equation on the brane |^|. 

The held equations induced on the brane are derived via an elegant geometric 
approach in [Q, leading to new terms that carry bulk ehects onto the brane: 


Gnu — + I^^Tnu + I^^Snv ~ £, 


■'flu ? 


(5.2) 


where = Sn/M^. The various energy scales are related to each other via 


K 


\ = 6'—, A = {A + 

Kj 


(5.3) 


and the high-energy regime is p ~ A. Bulk corrections to the Einstein equa¬ 
tions on the brane are of two forms: hrstly, the matter helds contribute local 
quadratic energy-momentum corrections via the tensor and secondly, there 
are nonlocal ehects from the free gravitational held in the bulk, transmitted via 
the projection of the bulk Weyl tensor. The matter corrections are given by 




flU 


= —T 

12 ' 


T 

a fiv 


Irji rpOi 
4 fia ^ u 


[3r„^T“'’ - (r„“)"] . 


The projection of the bulk Weyl tensor is 

£ab = CACBon^n ^, 


(5.4) 


(5.5) 


which is symmetric and traceless and without components orthogonal to the 
brane, so that Sab^^ = 0 and £ab —^ £iJLugA^gB'' as y —> 0. 

The Weyl tensor Cabcd represents the free, nonlocal gravitational held in the 
bulk, i.e., the part of the held that is not directly determined at each point by the 
energy-momentum tensor at that point. The local part of the bulk gravitational 
held is the Einstein tensor Gab^ which is determined locally via the bulk held 


equations (|5.1j). Thus £^u transmits nonlocal gravitational degrees of freedom 


from the bulk to the brane, including tidal (or Coulomb), gravito-magnetic and 
transverse traceless (gravitational wave) ehects 
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If is the 4-velocity comoving with matter (which we assume is a perfect 
fluid or minimally-coupled scalar held), we can decompose the nonlocal term as 

iu^iUu + -|- QuUf^ 1 (5-6) 

hv A 

where + u^Uy projects into the comoving rest-space. Here 

is an effective nonlocal energy density on the brane (which need not be positive), 
arising from the free gravitational held in the bulk. It carries Coulomb-type ehects 
from the bulk onto the brane. There is an ehective nonlocal anisotropic stress 

V^y = -\^X [h,^hy^ - 

on the brane, which carries Coulomb, gravito-magnetic and gravitational wave 
ehects of the free gravitational held in the bulk. The ehective nonlocal energy 
hux on the brane, 

Qf, = lK^Xhf,°‘Sa0U^ , 


carries Coulomb and gravito-magnetic ehects from the free gravitational held in 
the bulk. (Note that there is no energy hux in the bulk, and thus no transfer 
of energy between bulk and brane; this situation changes if bulk scalar helds are 
present |57|, PD|.) 


5.3 Local and nonlocal conservation equations 

The local and nonlocal bulk modihcations may be consolidated into an ehective 
total energy-momentum tensor: 

G,y = -Ag,y + KXu, (5-7) 

where 

T;Z' = T,y+^S,y-\£,y. (5.8) 

A tv 

The ehective total energy density, pressure, anisotropic stress and energy hux are 



/ p \ &J 

- ''(^ + 2a) + k-A ’ 

(5.9) 

ptot 


(5.10) 

tot 

6 

(5.11) 

tot 

V 

6 

(5.12) 
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The brane energy-momentum tensor separately satisfies the conservation equa¬ 
tions, = 0, giving 


p + 0(p + p)=O, (5.13) 

+ = 0 ) ( 5 - 14 ) 

where a dot denotes is the volume expansion rate of the 

congruence, is its 4-acceleration, and is the projected covariant 

spatial derivative. The Bianchi identities on the brane imply that the projected 
Weyl tensor obeys the constraint 


VI'S,,. = — . 


(5.16) 


This shows how nonlocal bulk effects are sourced by local bulk effects, which in¬ 
clude spatial gradients and time derivatives: evolution and inhomogeneity in the 
matter fields can generate nonlocal gravitational effects in the bulk, which back- 
react on the brane. The brane energy-momentum tensor and the consolidated 
effective energy-momentum tensor are both conserved separately. Conservation 


of gives, upon using Eqs. (^.9|) - (|5.14|) , propagation equations for the nonlocal 


energy density U and energy flux In linearized form, these are 


W + |0W + D'^Q^ = O, (5.16) 

+ + p)D/,p, (5.17) 

where H = a/a {= 40) is the Hubble rate in the background. The nonlocal 
tensor mode, which is the part of that satisfies = 07 ^ does not 

enter the nonlocal conservation equations. Furthermore, there is no evolution 
equation at all for reflecting the fact that in general the equations do not 
close on the brane, and one needs bulk equations to determine brane dynamics. 
There are bulk degrees of freedom whose impact on the brane cannot be predicted 
by brane observers, for example, incoming gravitational radiation from the bulk. 
The evolution of the nonlocal energy density and flux, which carry scalar and 
vector modes of the bulk gravitational field, is determined on the brane, while 
the evolution of the nonlocal anisotropic stress, which carries scalar, vector and 
tensor modes of the bulk field, is not. 

The generalized Raychaudhuri equation on the brane in linearized form is 


0 + 402 - + \k\p + 3p) - a 

=-4fi:2(2p + 3p)y -(5.18) 

A K A 

where the general relativistic case is recovered when the right-hand side is set to 
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zero. In the background, this gives 


H = 


K 


2 r 


P + 3p + ^(2p + 3p) 


-A 




Uo 


(5.19) 


where the solution for U follows from Eq. (|5.16|) , Oq is the initial scale factor 
and Uo = U{ao)- The hrst integral of this equation is the generalized Friedmann 
equation on the brane: 




P 


H - —p (1 + — )+-A- — + 


K 




(5.20) 


where A' = 0, ±1. Local bulk effects modify the background dynamics. In par¬ 
ticular, inflation at high energies (p ~ A) proceeds at a higher rate than the 
corresponding rate in general relativity. This introduces important changes to 
the dynamics of the early universe M, and accounts for an increase in 


the amplitude of scalar and tensor fluctuations at Hubble-crossing [^, ^ 
Using Eqs. ^.20| and |5.19| the condition for inflation becomes 


82 


1 / 2p -|- A 

< -77 — -W 

3 \ p -|- A 


(5.21) 


where w = p/p. As p/A —>• oo, we have w < —|, while the general relativity 
condition w < —| is recovered as p/A —>• 0. 

If Uo = 0, i.e., if the background bulk is conformally flat, then Eqs. (0) 
and (|5.10|) show that the effective equation of state index for the total energy- 
momentum tensor is 




w + (1 + 2w)p/2X 


P 


■tot 


1 + p/2A 


l + 2w, 


(5.22) 


where the last equality holds at very high energies (p S> A). Thus for slow-roll 
inflation, and w are both close to —1. The high-energy inflation condition 
to < — I is < — |- During high-energy reheating with tc 0 on average, we 
have ~ 1, so that the effective equation of state is stiff, while high-energy 
radiation-domination {w = |) has ~ |, i.e., an ultra-stiff effective equation 
of state. The effective sound speed at very high energies is also altered: 


where cf = p/p. 


Atot 




p 


(5.23) 
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5.4 Scalar perturbations on the brane 

When the vorticity of the fluid vanishes, scalar perturbations are covariantly (as 
well as gauge-invariantly and locally) characterized as the case when all per¬ 
turbed quantities are expressible as spatial gradients of scalars. In particular, 


the nonlocal perturbed bulk effects are described by and |]84 


Qm = , Pm- = [h/h/ - V] D„D^P . 

Note that the nonlocal traceless mode has spatial divergence ^4 


(5.24) 


= |D2(D^P) . 

The bulk gravitational field affects scalar perturbations via scalar Coulomb modes, 
given by the spatial gradients of the ‘potentials’ U, Q and V. 

For adiabatic matter perturbations the 4-acceleration is 




p{l + w) 


DmP- 


(5.25) 


The gradients 




(5.26) 


describe inhomogeneities in the matter and expansion and the dimensionless 
gradients describing inhomogeneity in the nonlocal quantities are 


, Q^ = -B^Q, P, = —D,V . 
p p ap 


(5.27) 


The spatial gradient of the conservation equations ( |5.13|) , ( p.l6|) and (|5.17|) , 
and the generalized Raychaudhuri equation (|5.18|) , leads to a system of equations 
for these gradient quantities |^. The gradients define scalars via their comoving 
divergences: 


F = , with F = A,Z,U,Q,P, (5.28) 

where A is a covariant analog of the Bardeen density perturbation Cm (see Eq. |2.1CI|) 
Then the system of equations governing scalar perturbations on the brane follows 
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from the gradient system given in as 


A = 3wHA — (1 + w)Z , 


Z = -2HZ - 


1 + w 


D^A - 


6p 


U 


Wp 


P 


l + (4 + 3w)^ 


2c? \ 6U 


A \1 + w J K‘^Xp 


A, 


(5.29) 


(5.30) 


il = {3w-l)HU- 


4c2 


1 + ta 


HA 


A^Z-aD^Q, 


g = (1 - 3w)HQ - U - laD^P 

3a 


Sd 


6a L\l + ta/ p 


U 


— K^p{l + w) 


A. 


(5.31) 


(5.32) 


In general relativity, only the hrst two equations apply, with A“^ set to zero in 
Eq. ( p.30| ). In this case we can decouple the density perturbations via a second- 
order equation for A, whose independent solutions are adiabatic growing and 
decaying modes. Local bulk effects modify the background dynamics, while non¬ 
local bulk effects introduce new fluctuations. This leads to fundamental changes 
to the simple general relativity picture. There is no equation for P, so that in 
general, scalar perturbations on the brane cannot be predicted by brane observers 
without additional information from the unobservable bulk. Thus in general, one 
must solve also the scalar perturbations in the bulk in order to determine the 
perturbation evolution on the brane. 

However, there is a crucially important exception to this, arising from the 
fact that P only occurs in Eqs. (|5.29|) - (|5.32|) via the Laplacian term D^P. We 
can use the shear propagation equation on the brane [Q to provide an order-of- 
magnitude comparison of P with A: 


a 


fJ^lP 








_p 


where is the electric Weyl tensor on the brane. This equation, together with 
Eqs. (^.24|) and (|5.25|) shows that 




ID^PI 


1 

P' 


dVI, 


and then Eqs. (|5.26|) - (|5.28|) imply 


IPI 


K^X 

a?p 


|A|. 
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Then 


4 

K P\ 




a?H^ a ' 


on using the high-energy Friedmann equation ~ /\). Thus, for k -C aH, 

i.e. on large scales, well beyond the Hubble horizon, we can neglect the D^P 
term in Eq. (|5.32|) relative to the A term. Thus on large scales, the system closes 
on the brane, and brane observers can predict scalar perturbations from initial 
conditions intrinsic to the brane, without the need to solve the bulk perturbation 
equations. Note that the term in Eq. (|5.31|) may also be neglected relative 
to the U term. This follows from the shear constraint equation 




which gives, on taking the divergence. 


IQI 


ap an 


where the last relation follows from Eq. (|5.3CI|) . Thus 


\aD^Q\ 


a2p2 


{H\U\) 


and so we may neglect the D^Q term in Eq. |5.31| relative to the U term. The 
system Eqs. (|5.29|) - (|5.32|) then reduces to 3 coupled equations in A, Z and U, plus 
a decoupled equation for Q, which determines Q once the other 3 quantities are 
solved for. Thus there are in general 3 modes of large-scale density perturbations: 
a non-adiabatic mode is introduced by bulk effects. This mode is carried by 
fluctuations U in the nonlocal energy density U, which are present even if U 
vanishes in the background. The fluctuations Q and P in the nonlocal energy 
flux and anisotropic stress do not affect the density perturbations on very large 
scales. 

In qualitative terms, we can interpret these general results as follows. Bulk 
effects lead to an effective total energy-momentum tensor that is non-adiabatic. 
From Eqs. ( p.9|) and (|5.10| ), we hnd a measure of the total effective non-adiabatic 
pressure perturbation, which is the covariant analog of 


aDf,p 


tot 


X 


3 


p + p 


6Hp 




p + X 


^4ptoty 

{4UA^-3{p + p)U^} . 


(5.33) 


In addition, Eqs. (|5.11|) and (|5.12|) show that non-adiabatic stresses and fluxes are 
also present, due to nonlocal bulk effects. From this viewpoint, it is not surprising 
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that non-adiabatic modes arise in the density perturbations. An analysis of the 
non-adiabatic large scale mode in terms of the curvature perturbation is given in 


S5 


Note that an alternative interpretation is also possible. We can regard the 
nonlocal bulk effects as constituting a radiative “Weyl” fluid, with energy-momentum 
tensor This Weyl fluid has non-adiabatic stresses and is in mo¬ 

tion relative to the matter fluid, with relative velocity parallel to Q^. Although 
the matter fluid obeys energy-momentum conservation, the Weyl fluid does not; 
the momentum balance equation ( ^.171) shows that density inhomogeneity in the 
matter fluid sources a momentum transfer to the Weyl fluid. 

Whatever intuitive interpretation we adopt, the concrete result is that large- 
scale density perturbations on the brane may be determined without knowledge 
of the bulk, and acquire a non-adiabatic mode due to effects from the free gravi¬ 
tational field in the bulk. 


5.5 Large-scale density perturbations 

We rewrite the coupled system for large-scale perturbations by introducing two 
useful new quantities. We dehne, following |^ , 


= K^pa^A , 


(5.34) 


which is a covariant analog of the Bardeen metric potential (see Chapter 
and the covariant local curvature perturbation 


C = , Cy = a^DyR^ , 


(5.35) 


where is the Ricci curvature of the surfaces orthogonal to u^. (Note that 
these surfaces are in general shearing, and non-uniform in p, Q, U and R^.) 

Then the coupled system for density perturbations, Eqs. (|5.29|) -( p.31|) , can 
be rewritten on large scales as 


$ = 


-H 


(1 + w)kV / p 

2H^ V A 


$ 


C = - 
U = H 


(1 -|- w)h?p 
TR 
72clHU 


(1 -|- w)XiRp 


C - 


3(1 -I- w)a^p^ 

w 


u, 




3w-l- 


4U 


2U 

3aXHp 


i + ^ + 




QclH^ 


X {l + w)iRp 




(5.36) 

(5.37) 


(5.38) 


51 

































The general relativistic case is recovered when we set A W and U to zero; 
in this case, Eq. (p.38|) falls away, and Eq. ( ^.371) reduces to 


^ = ^0, a = o. 


(5.39) 


which expresses conservation of the covariant curvature perturbation along each 
fundamental world-line. The value of Co will in general vary from world-line to 
word-line, so that its conservation is local, and is not an indicator of purely 
adiabatic perturbations. (In general relativity, C = 0 on large scales for a flat 
background even when there are non-adiabatic perturbations [^.) Bulk effects 
destroy the local conservation of C in general, by Eq. (|5.37|) . 

However, there is an important special case when local conservation is re¬ 
gained: when the nonlocal energy density U vanishes in the background. This 
does not mean that fluctuations in the nonlocal energy density are zero, i.e., we 
still have 17 7 ^ 0 in general. It can be argued that vanishing U in the background 
is more natural, if one believes that the bulk background should be conformally 
flat, and thus strictly anti-de Sitter. (Quantum effects may nucleate a black 
hole in the bulk |^, in which case the Schwarzschild-anti de Sitter bulk, with 
Uo proportional to the black hole mass, would be a natural background.) From 
now on, we will assume a conformally flat background bulk and a spatially flat 
brane background; thus Uo = 0 = K in the background generalized Friedmann 
equation (|5.20|) . Equation (|5.33|) shows that the non-adiabatic total pressure per¬ 
turbation is then proportional to {p/X)U, which will be enhanced at high energies 
and suppressed at low energies. 

When U = 0 in the background, Eq. (|5.39|) holds, and Eq. (|5.38|) gives 


U = Uof , ilo = 0 , f = exp / {3w — 1) d In a . 


(5.40) 


This shows that U rapidly redshifts away during inflation, so that non-adiabatic 
effects from nonlocal bulk influence are small. By contrast, the modifications to 
the background dynamics from local bulk effects are strong during inflation at 
high energy. 

The key equation (|5.36|) becomes 


dd) 

dN 


1 + 


(1 -I- w)k^p 




$ = 


(1 -I- w)iU‘p 

Co - 

3(1 w)alp^ 

[ 4772 J 


[ \ 


-,2N 


fUo, 


(5.41) 


where N = ln(a/ao) is the number of e-folds. We have thus reduced the coupled 
system to one simple inhomogeneous linear equation, which may be integrated 
along the fundamental world-lines. Along each world-line, the constancy of Co 
and Uo allows us to track the change in $ as tc changes, from inflationary behavior 
through to radiation- and matter-domination. 
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We can perform a qualitative analysis of the evolution of as follows. Substi¬ 
tuting from the generalized Friedman Eq. ( ^-201) into Eq. ( b.41|) and setting 
Uq, K and A to zero gives 


dN 


$ 


+ (l + I (1 +^) (l + a) + I a) ) 

= f (1 + ^) C'o (l + I y) 

-9 (1 + w) (l + H) 


(5.42) 


-1 


For high-energy inflation p S> A and so we can approximate Eq. (|5.42|) by 

(l + u;)ao2e^(^+3“')t/o 


(i<F , ^ ^ 3 (1 -|- tc) Cq\ 
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(5.43) 


As it is slow-roll inflation, we can assume w and p are constant when solving this 
equation. This gives 


<h (1 -I- tc) - 


ACn 


2 (4 -I- 3tc) p 


- 18 


(5 -|- 6 ta) 


^ + g(-4-34«)iv^^ ^5 


44) 


where Ci is an integration constant. Using w 
tially decaying terms this gives 

high-energy inflation: <F 


T and dropping the exponen- 


-(1 w)Co — . 

2 p 


(5.45) 


For the general relativity limit we use A S> p and then following the same proce¬ 
dure as for the high energy case we get 


~ + w)Co 


(5.46) 


In general relativity, <F remains constant on large scales during slow-roll inflation, 
independent of the form of the inflaton potential. In the brane-world, <I> is slowly 
increasing during high-energy slow-roll inflation, since <I> ~ p~^ and p is slowly 
decreasing. This qualitative analysis is conhrmed by the numerical integration of 
a simple phenomenological model shown in Fig. 1. We have modelled a smooth 
transition from inflation to radiation by tc = 4[(2 — a) tanh(W — 50) — (1 — a)], 
where a is a small positive parameter (chosen as a = 0.1 in the plot). For more 
realistic models, i.e., where U((p) is specihed, the evolution of <I> may be more 
complicated than shown in Fig. 1. For Po/A S> 1, Eq. (p.21| ) shows that inflation 
ends at W = 50 — 21n[(l — 2a)/3] ~ 47.4, and at W = 50 in general relativity. 
Only the lowest curve still has p/A 3> 1 at the start of radiation-domination [N 
greater than about 53), and one can see that <h is still growing, as conhrmed by 
Eq. (1^) . 
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During reheating, in periods when w is approximately constant on average 
(for example, ta 0 for D = Eqs. ( ^.131 ) and (|5.41|) imply 


high-energy w ~ constant reheating: 

® ^ + const. (5.47) 

2{7 + 6w) po 

Thus high-energy w ~ constant reheating on the brane produces amplihcation 
of <h, unlike general relativity, where $ remains constant on large scales during 
w ~ constant reheating: 


$ 


gr 


3(1 + ^) 

2(5 + 3m;) ' 


(5.48) 


In the radiation era, the energy density redshifts rapidly, so that p quickly 
falls below the brane tension A. If the energy density at the end of reheating is 
high enough, then at the start of radiation-domination we have p 3> A, and we 
hnd that <h is amplihed during high-energy radiation domination: 

2 A 

high-energy radiation: $ - Coe^^ -\- const. (5.49) 

9po 

At low energies on the brane, or in general relativity, we hnd that $ is constant: 


low-energy radiation: 


$ $gr 



(5.50) 


This qualitative result is conhrmed in Fig. 1. After the radiation era, the energy 
scale has fallen well below the brane tension, so that in the matter era, we recover 
the general relativity result: 


matter era: 


4*. « IqC"- 


(5.51) 


In general relativity, the constancy of $ during slow-roll inhation and radiation- 
and matter-domination allows one to estimate the amplihcation in $. CMB large- 
angle anisotropies as measured by COBE place limits on the amplihed <h, and 
this in turn places constraints on the inhationary potential, since the potential 
determines the initial value of <h. This simple picture is complicated by high- 
energy and non-local ehects in the brane-world. Further discussion of large-angle 


CMB anisotropies on the brane is given in |85 . 


5.6 Conclusions 


Using the 1-1-3 covariant formalism developed in [Q, we have analyzed the evo¬ 


lution of large-scale density perturbations on the brane. Density inhomogeneity 
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Figure 5.1: The evolution of $ along a fundamental world-line for a mode that 
is well beyond the Hubble horizon at iV = 0, about 50 e-folds before inflation 
ends, and remains super-Hubble through the radiation era. Labels on the curves 
indicate the value of po/ A, so that the general relativistic solution is the dashed 
curve (po/A = 0). 
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on the brane generates Weyl curvature in the bulk, which in turn backreacts on 
the brane, in the form of a nonlocal energy-momentum tensor. Fluctuations in 
the nonlocal energy density induce a non-adiabatic mode in large-scale density 
perturbations. The fluctuations in the nonlocal energy flux are decoupled from 
the density perturbations, while the nonlocal anisotropic stress plays no role on 
large scales. This latter feature is what closes the system of brane density per¬ 
turbation equations, allowing brane observers to evaluate the perturbations on 
the brane without solving for the bulk perturbations. 

We showed that the local and nonlocal bulk effects arising during high-energy 
inflation, and any high-energy start to the radiation era, modify the simple pic¬ 
ture of general relativity. The local covariant version of the metric perturbation, 
i.e., <F, is no longer constant on large scales during these regimes. We gave a 
rough estimate for slow-roll inflation in Eq. (^.45|) . Numerical integration will 
be required for the more complicated case, not investigated here, when the back¬ 
ground bulk is not conformally flat, i.e., when the nonlocal energy density U 
does not vanish in the background. In this case, the coupled system of equa¬ 
tions ( b-36|) -( p.38| ) can no longer be reduced to one equation, since C and U are 
no longer locally conserved. 

The formalism we have used is restricted to large scales. When a mode ap¬ 
proaches the Hubble radius, the gradient terms can no longer be neglected, and 
the presence of these terms means that the system of equations no longer closes on 
the brane. A fuller investigation requires a formalism that can handle all scales, 
and which necessarily involves the evolution of perturbations in the bulk. A co¬ 
variant formalism for bulk cosmological perturbations has not been developed, 
but a metric-based formalism has been developed |5B, 0- The equations of this 
formalism are very complicated, and considerable work remains to be done before 
smaller scale structure can be predicted and compared with observations of the 
acoustic peaks in the CMB anisotropies. Our results provide a useful initial step 
for further developments by showing what happens on very large scales. 
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Chapter 6 


The effect of non-linear 
inhomogeneities on inflation 

6.1 Introduction 

In this Chapter we discuss in more detail the initial conditions required for in¬ 
flation. In particular we examine what scales of non-linear perturbations can be 
smoothed out by inflation. 

Inflation is the foremost idea for explaining the large scale homogeneity and 
near flatness of the universe, which are the two main unexplained observational 
features in the standard hot big-bang model. The large scale homogeneity or 
horizon problem amounts to the fact that under hot big-bang evolution, due 
to a decelerating scale factor d < 0, sufficiently separated regions of the present 
day observable universe would never have been in causal contact. Nevertheless, it 
appears the universe we observe looks very much the same, and in particular very 
smooth, in all directions. This fact is best seen in the CMB which has temperature 
fluctuations of only one part in 10^ when measured from any direction in the sky 

i- 

The inflation solution to this horizon problem is to picture the universe during 
an early epoch to undergo an accelerated expansion, d > 0. Such expansion can 
take an initially small causally connected patch and enlarge it to a size that 
comfortably encompasses our present day observed universe, thereby solving the 
horizon problem. To realize inflation, the equation of state within the inflationary 
patch must be of a very special form, possessing negative pressure. The potential 
energy of a scalar field has an equation of state that satisfies this requirement. 
This fact has been a key link towards a dynamical realization of inflation and 
thereby has further motivated the inflation solution. 

As inflation is meant to solve the horizon problem, it is important that it 
does not require acausal initial conditions. In particular, the picture of inflation 
considered here is for the universe to emerge from an initial singularity and then 


57 


enter into a hot big-bang radiation dominated evolution. At some time ti after 
the initial singularity, the conditions appropriate for inflation should occur within 
a small patch that is contained within the causal horizon at that time. Chaotic 
inflation ^ does not fall into this picture as there inflation is thought to 
start at the Planck epoch with homogeneity assumed on the Planck scale. To 
realize the picture of a “local” inflation, two requirements must be satished. First 
a physically sensible embedding must be demonstrated of the inflating patch 
immersed within a non-inflating background |^, |100|| . By embedding we 

mean matching the inflationary space time with the background space time at 
the boundary of the patch. Second, it must be shown that for t > ti, the patch 
is dynamically stable to sustain inflation ||101| , |102| , |103|, |104|, p.05| , p.06| , |107| . For 
example, large fluctuations, which conceivably could enter the inflating patch at 
a maximum rate limited only by causality, should not destroy the inflationary 
conditions within the patch. 

Recently, a convenient methodology has been developed in |^, |100|| for analyz¬ 
ing the embedding problem, based on the null Raychaudhuri equation ||1U8| , 

In Sec. |6.2| the flat spacetime (hi = 1) formulation in 
arbitrary spacetime curvature. Then, an alternative solution from 


is generalized to 
to the 


embedding problem will be identihed, which is especially attractive for an in¬ 
flating patch with an open geometry. In Sec. |6.3| , initial conditions for scalar 
held dynamics are presented, which are consistent with causality and our embed¬ 
ding solution and which evolve into successful supercooled or warm inhationary 
regimes. In Sec. we combine the results in Sec. and Sec. to evalu¬ 
ate the ehect of inhomogeneities on the embedding problem. Finally, Sec. 
presents our conclusions. 


6.2 Embedding conditions 


The embedding of an inhationary patch within a background space time gener¬ 
ally is regarded as unacceptable if there are any negative energy regions. This 
requirement often is referred to as the weak energy condition. The null Ray¬ 
chaudhuri equation is a useful diagnostic tool for determining the validity of this 
condition. This equation determines the evolution of the divergence 6 = V“Aq, 
of the null ray vector in a spacetime with an arbitrary, and not necessarily 
homogeneous, energy density distribution. In order for the weak energy condition 
to be valid, this equation implies that for a null geodesic the condition |^, |109 


d9 


( 6 . 1 ) 


must be satished, where s is an affine parameter along the null geodesic. 

For application of Eq. (|6.1|) to the inhation embedding problem, the concept 
of anti-trapped and normal regions should be understood. Consider a sphere 













































centered on a conioving observer. If the space time were not expanding, photons 
emitted from the surface of the sphere, radially towards the observer, would 
converge at the observer. However, the expansion tends to work against the 
bundle of rays converging to a point. If the expansion is rapid enough, the bundle 
of rays will have diverging trajectories and then the spherical surface from which 
the rays originated is said to be an anti-trapped surface [|I09|| . The spherical 
surface with a radius Xmas is known as the minimally anti-trapped surface (MAS) 
if any sphere with a larger radius is anti-trapped. For inwardly directed null rays, 
the divergence, 6, will be positive in anti-trapped regions of space. On the other 
hand, if 6 is negative for inwardly directed null rays and positive for outwardly 
directed null rays, the region is called normal. 

These dehnitions are convenient when analyzing the weak energy condition 
based on Eq. (O). For example, one can immediately conclude [^| that if 
an outer normal region bounds an inner anti-trapped region in a not necessarily 
homogeneous, spherically symmetric space time, the weak energy condition would 
be violated, thus implying negative energy is required at the boundary of such a 
conhguration. 

For the inflation problem, the inner region is modeled as the putative infla¬ 
tion patch (INF) and it is immersed within a outer background (BG) expanding 
spacetime region. As a first approximation, both the INF and BG regions are 
characterized by FRW metrics of the form 


ds^ = —dr -I- 


dr‘^ 


l-Kr^ 


+ r^dVL^ 


( 6 . 2 ) 


where K = —1 for an open geometry, 0 for a flat geometry and 1 for a closed 
geometry. The effect of inhomogeneities in the BG will be considered in Sec. 
6.4 . The INF patch should be homogeneous. In what follows, 9 and Xmas will 
be computed for a spacetime characterized by the metric Eq. ( |6.2|) . The results 
below are applicable for both the INF and BG regions, given that the appropriate 
parameters are used in Eq. ( |6.2|) for the two different spacetime regions. 

Proceeding with this calculation, the determinant of the metric Eq. (|6.2|) is 


9= |det(g„„)| = oV(l - AV“) 


(6.3) 


where is the polar angle. An incoming light ray (null vector) is given by 


AT” = i [dj - a-‘(l - A'r^)‘/^(5f] , 


(6.4) 


with the divergence of the null rays being 


1 


9 = 

V9 


(6.5) 
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Substituting equations (|6.3| ) and ( |6.4|) into ( |6.5|) gives 


9 = 


H- 




ar 


( 6 . 6 ) 


with H = a/a the Hubble parameter. The above recovers equation (10) of [TO 
for = 0. 

The physical distance is given by 


X = a 


,2^-l/2 

(6.7) 

for K = -l 

(6.8) 


The MAS comoving distance, must give a zero divergence and so from Eq. 
with, for example, K = —1 


H 


("X 7. 2\l/2 

yj. -r I xnas J 


= 0 , 


mas 

from which we get, using equation (|6.^) , 

X mas = aarcsinh [(H^a^ — 1)”^^^] . 
From the Friedman equation: 

K 


H^a^ = 


n 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


Substituting (|6.11|) into (|6.10|) gives the solution for Xmas as a function of H < 1. 
A similar equation can be derived for the closed case (1 < H < 2) where an upper 
limit is needed on H to ensure that 1/H is smaller than the radius of the Universe. 
Therefore, the general MAS size is given by 







( 6 . 12 ) 


A plot of equation (|6.12|) is given in Fig. |6]^. As can be seen, the size of the 
MAS becomes arbitrarily large relative to the Hubble radius as H ^ 0. 

The results for the Xmas derived above are valid for either the INF or BG 
region. In [Q, the case H = 1 was treated, and our calculation above is consistent 
with them in this limit. Before proceeding with our analysis for general H, let us 
review the results of for H = 1. The background spacetime in their work is 
pictured, as for us, as evolving in a hot big bang regime. For H = 1 and standard 
forms of matter, e.g. radiation, non-relativitic matter, etc., the Hubble horizon 
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Figure 6.1: A plot of MAS size as function of the density. 


s 



R 



Figure 6.2: A Penrose diagram for local inflation (adapted from ). The arrow 
(ab) denotes a radially directed null-geodesic going from the normal background 
space time to the anti-trapped space time inside the patch. Shading represents 
anti-trapped regions. 
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sets the the scale on which causal processes can take place. For dehniteness, if 
we take the background causal horizon size as Xp® = 1/Hbg, then by Eq. (|6.12|) 
for = 1, Xp® = xj^as- Fig- |6.2| shows a Penrose diagram for an embedding that 


violates the weak energy condition. Here, at the beginning of inflation, Xp is 
represented by line OQ and x™^g is represented by line OP. An inflating patch is 
pictured to develop within some region inside the background (polygon OQRS in 
Fig. |6]^). In the inflationary patch, they assume a different Hubble parameter 
Hinf and assume the size of this patch must be x™^ > l/Hinf. This assumption can 
be justihed by requiring the patch to be stable against perturbations and will be 
discussed further in the next section. By Eq. (|6.12|) it also implies x™^g = l/Hinf 
for O = 1. If the patch is to be set up by causal processes then we would expect 
it to be smaller than the background Hubble horizon, x™^ < 1/Hbg. These 
conditions combined imply x((fns = 1/Hbg > 2 :™^ > x^^s- As such the region 


is normal with respect to the BG-region and the region 
is anti-trapped with respect to the inflating patch. Thus an 


between xjjf^s ^md Xp 
between x|f^ and x“L 
in-going null ray will go from a negative divergence region to a positive divergence 
region. Based on Eq. (|6.1|) this leads to a violation of the weak energy condition. 
Due to this fact, in they conclude that the only way to avoid this violation is 
to have x™^ > 1/Hbg which may be impossible or at least very difficult to come 
about through causal processes. 

Although this analysis assumed a homogeneous background, the general result 
of [^] was that if the inflationary patch contained a MAS then it must be larger 
than the background MAS. One of the main aims of this article is to further 
examine the implications this has for a causal embedding of an inflationary patch. 

Our Erst observation is that there is an embedding consistent with both 
causality and the weak energy condition, in particular 


/-y> 

5 ‘^mas ^ 


inf 


(6.13) 


Since x™^ in Eq. ( |6.13| ) is smaller than Xp®, it implies consistency with causality. A 
region which is smaller than its MAS will have no anti-trapped surfaces. Thus in 
our proposed embedding Eq. (|6.13|) , the INF-BG boundary is between two normal 
regions. As such, provided 6*™^]boundary is sufficiently smaller than 6*’’®]boundary, the 
weak energy condition is satisfied. One interesting feature of the embedding Eq. 
(|6.13|) is that for Ojnf < 1 it is acceptable for x™^ > l/Hj, 


Mnf- 


As inflation proceeds, Oinf —>■ 1 and so x™^g —l/Hjnf. One of the main 
features of inflation is that modes of a matter perturbation become larger than 
the Hubble radius as the Universe expands. It follows that eventually x“^ > 


inf 


l/Hjnf, X 

condition provided that x™^ > xlij®, at that time 


Fig. 


will need to occur. This will not cause violation of the weak energy 

inf 
P 

3| shows a Penrose diagram of the proposed embedding. At point T 


in the hgure, x™^ = x|(fa^g and at point P, 


X 


inf 


= X 


inf 


As can be seen the anti- 


p mas ^ 

trapped region develops in the patch only after the patch has become greater 
than the background MAS. 
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Figure 6.3: A Penrose diagram for local inflation with an embedding that does 
not violate the weak energy condition and is not acausah 


There is another embedding which does not violate the weak energy condition, 
in this case for ffbg > 1, but it proves not to be as useful, 

rbg 


™inf ™bg ^ ™inf 
*^mas5‘^mas ^ *^p 


< l/Hbg. 


(6.14) 


This solution is restricted by the upper limit on flbg that gives a lower limit (from 
Eq. |6.12|) of 0.8/Hbg for = 2. As the causal scale can be only of the 

order of l/Hbg, such a large inflationary patch may still be acausah 

Additional information about embedding constraints can be obtained through 
the Israel junction conditions [|110|| . In [|111|, [112|| conditions were derived for the 


embedding of one FRW spacetime within another. The Israel junction conditions 
require the extrinsic curvature on the background side of the boundary to be 
smaller than the extrinsic curvature on the patch side of the boundary if the 
boundary surface energy density is to be positive. It can then be shown ||111|, |112| 


that if the energy density of the background pbg is greater than that of the patch 
Pinf, the Israel junction conditions can be satished for an inflationary patch of 
arbitrary size and geometry with any background geometry. They also show 
for pbg < Pinf, an embedding consistent with the Israel junction conditions and 
< l/Hbg requires the background geometry to be closed with the geometry 
of the inflationary patch arbitrary. 
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6.3 Dynamic Conditions 

An immediate concern with the proposed embedding Eq. (|6.13|) regards the size 
of the inflationary patch. As far as the embedding conditions are concerned, 
the patch size simply must lie below the x^as line in Fig. 1. In particular, for 
riinf > 1 the inflationary patch size must be < 1/Hinf. However, for Hinf < 1, the 
inflationary patch size can be > 1/Hinf. The primary question which remains is 
what the minimal initial inflation patch size can be in order to be dynamically 
stable for inflation to commence. For example, if the initial patch is too small, 
then large fluctuations initially outside the patch could enter inside at a rate 
sufficiently fast to destroy the inflationary conditions. In this section, the dynamic 
conditions necessary for inflation will be examined for both the scalar held and 
a background radiation component. 


6.3.1 Scalar field 


For scalar held driven inhation, the general dynamical conditions necessary for 
a hnite spatial patch to initiate and sustain inhation have been considered in 
101| , |102| , |103i |104| , |105| , |106|| and a comprehensive review has been given in [p.07 


Here the requirements addressed in these works, primarily the review |p.07|| , will be 
systematically examined with emphasis on determining their implications for the 
minimal initial inhationary patch size. Our considerations of scalar held dynamics 
will generalize those in the above stated works, which focused on supercooled 
inhationary dynamics [nB, ni, in that here warm inhation dynamics [ 115 
also will be treated. 

The classical evolution equation for the scalar inhaton held has the general 
form 


0(x, t) + [3H + r]0(x, t) 

- a~ V V(x, t) + = 0, (6.15) 

where T is a dissipative coefficient, which represents the ehective interaction of the 
inhaton with other helds. Here and below 0(x, t) is a classical real number which 
contains the ensemble average of the quantum and/or thermal huctuations. In 
supercooled inhation, it is assumed the inhaton is isolated, in which case T = 0. 
Thus no radiation is produced during inhation and the universe inhates in a 
supercooled state. On the other hand, in warm inhation the inhaton interacts 
with other helds, and the Tcp term is the simplest representation of their ehect 
on the evolution of the inhaton, with generally T > H. Due to these interactions, 
radiation is dissipated from the inhaton system into the universe throughout the 
inhation period. In general, T can vary as a function of the inhaton held mode 
to which it is associated, but here T is treated as a constant. 
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There are two types of potential V (0) in Eq. (|6.15|) that typically are treated 
in inflationary cosmology, a purely concave potential of the generic form V ~ 
and a double well potential such as E ~ (0^ — m^)^. Our discussion 
to follow will focus of the former type of potential and at the end we will comment 
on the case of the double well potential. Furthermore, for most of our discussion 
it will be adequate to consider the simplest case of a quadratic potential, V = 
In this case, going to Fourier space where we put the universe in a box 

t) = (6.16) 

k 


Eq. (|6.15|) becomes, 

0k + [3H + r]0k(^) + a~‘^k‘^(j)k{t) + m^0k(^) = 0, 


(6.17) 


where = |kp with k the comoving wave-vector and k/a(t) the corresponding 
physical wave-vector at time t. 

A necessary condition for inflation is that the zero mode 0k=O; must have a 
sufficiently large and long-sustained amplitude so that the potential energy V = 
dominates the equation of state of the universe, thereby driving inflation. 
This leads to the familiar slow-roll conditions, which require the curvature of the 
potential to be sufficiently flat so that Eq. ( 6.17|) for the zero mode becomes first 
order in time 


dV{(p)/d(po 

3H + r ■ 


(6.18) 


The dynamic initial condition problem is that the above requirements should 
not be too special and in particular should not violate causality. The most acute 
initial condition problem discussed in |1U7| and related works is that generally the 
initial inflaton field configuration will be very inhomogeneous, thus considerable 
contribution of gradient energy ((V0)^) should be present. In its own right, the 
gradient term for comoving mode k has energy density pv = (k^/2a^)0k 
equation of state pv = —Pv/3. From the scale factor equation 


d AttG , , , , 

- = -^P + 3p, 6.19 

a 6 

the effect of the gradient energy vanishes on the right hand side. As such, any 
additional contribution from vacuum energy still should drive inflation. However, 
excited modes with k/a > 3H -I- T not only will possess gradient energy, but 
based on Eq. (|6.17|) are under-damped. As such, they also have a kinetic energy 
contribution ~ 0k) which has an equation of state = p^. If the kinetic 
energy components of these modes dominates the energy density in the universe, 
then from Eq. (p.l9|) inflation will cease to occur. 
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For the case of supercooled inflation, since F = 0, the simplest way to avoid 
this problem is to require that modes with fc/a ~ H initially should not be excited. 
In other words, the inflaton field initially should be smooth up to physical scales 
larger than ~ 1/H. However, under general conditions, the causal size of the 
pre-inflation patch also will be of order the Hubble radius 1/H. As such this 
homogeneity requirement on the initial inflaton held impinges on being acausal, 
since this condition essentially requires initially smooth conditions up to the 
causal scale. 

To treat excited modes with k/a > H in the supercooled inhation case, 
the initial condition dynamics are much more complicated. A simple analytical 
method applied to this situation is the effective density approximation presented 
105|| and reviewed in Sec. 7.2 of ||107|| . In this approach, the inhomo- 


m 


101 


geneities of the inhaton are treated in determining its evolution, but the effect 
of these inhomogeneities on the metric is only treated in the Friedmann equa¬ 
tion through homogeneous terms that represent the effective gradient and kinetic 
energy densities as 


a 




-0Q -I- -m^ 
2 


+ Pr 


+ 


K 


( 6 . 20 ) 


For initial held conhgurations with sufficiently excited fc/a > H modes so that 
their gradient energies dominate the equation of state of the universe, the ehective 
density approximation indicates the Friedmann and inhaton evolution equations 
are highly coupled. In particular, the Hubble parameter will be dominated by the 
gradient energy term and this will act back on the inhaton evolution through the 
3H0 term. Furthermore, the initially large gradient terms also can in turn induce 
large kinetic energy in the modes. The outcome found in (see also 

p.03|| ) for this case is that the universe expansion is non-inhationary. On the other 


hand, for the excited fc/a > H modes with smaller amplitudes, their evolution 
will be oscillatory and once again the universe expansion is non-inhationary. For 
either of these two possibilities, this approximation method hnds that after an 
initial period of detaining the universe in a non-inhationary regime, the ehect 
of these higher modes becomes negligible. At this point the remaining vacuum 
energy dominates and eventually inhation proceeds. 

|104| , [I05|| which exactly treat the ehects on the metric 


Numerical simulations 


due to the inhomogeneities in the inhaton held, in fact, do not support this hnal 
conclusion of the ehective density approximation. Instead, the simulations hnd 
that once modes with fc/a > H are sufficiently excited, inhation generally does 
not occur or is highly suppressed. Although the hnal conclusion of the ehective 
density approximation was incorrect, it at least indicated that there is nontrivial 
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interplay between the inflaton and metric evolution equations once significant 
short distance inhomogeneities are present in the infiaton field. 

In contrast to the supercooled case, once F > H, the effective density approx¬ 
imation indicates that for modes with F > fc/a > H, there is almost no feedback 
between the inflaton and Friedmann equations. The evolution equation for these 
inflaton modes is 


TMt) 





( 6 . 21 ) 


and in particular is independent of H. This in turn implies all these modes are 
over-damped. As such, they only contribute gradient energy to the equation of 
state of the universe, which alone is ineffective in preventing inflation. Thus up 
to the predictions of the effective density approximation, in contrast to the super¬ 
cooled case, in this warm inflation case modes with k/a <r show little coupling 
between the inflaton and Friedmann equations. It is evident that dynamics with 
a F0 damping term should have qualitative differences for the initial condition 
problem compared to the case with F = 0. In particular, modes of physical 
wavelength smaller that the Hubble radius 1/H but larger than 1/F could still 
be substantially excited without destroying entrance into inflation. This implies 
inflationary patch sizes smaller than the Hubble radius may be dynamically sta¬ 
ble in sustaining inflation and based on the results of Sec. |6.2| , they also provide 
consistent spacetime embeddings, especially for Hjnf < 1. Furthermore, studies of 
warm inflation ||115| , ^,|n3, including the first principles quantum field theory 
model ||118|| , generally find that to obtain adequate inflationary e-folds, ~ 60, 
it requires F Nem^/R with m ~ H so that F ~ A^"eH and thus 1/F -C 1/H. 
As such, for warm inflationary conditions this simple analysis suggests that the 
smoothness requirement on the initial inflationary patch is at scales much smaller 
than the Hubble radius 1/H, which therefore imposes no violation of causality. 

These considerations can be extended to the more complicated situation where 
the infiaton evolution equation has nonlinear terms, such as the interaction 
term. For supercooled inflation, the effect of such interactions has been consid¬ 
ered from the perspective of mode mixing for some special cases in 


102, 103 


and through computer simulations of scalar fields in one spatial dimension in 
1 1105|| . Up to the scope of these works, their analysis concludes that nonlin¬ 


ear interactions do not present any additional complications to the initial con¬ 
dition problem. For the warm inflation case, this conclusion can be stated in 
more general terms. In particular, the dissipative coefficient F completely damps 
the evolution, thus suppressing mode mixing, of all modes with k/a < F. As 
such, initially excited modes with k/a <T will evolve independently within the 
time scales relevant to the initial period when inflation begins. The only effect of 
excited modes with /c/a < F is to contribute gradient energy to the equation of 
state of the universe, and this alone is ineffective in suppressing inflation. There¬ 
fore, provided the inflaton field also initially maintains some vacuum energy, no 
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aspect of the inflaton’s dynamics in the early period acts to circnmvent entrance 
into the inflation phase. 

Although the above discussion addressed the main impediment in scalar held 
dynamics that can prevent inhation, there are a few smaller concerns worth men¬ 
tioning here. One detail treated in |[107]| is inhaton initial conditions with large 
kinetic energy. For F = 0, h is shown in |p.07|| that this problem is self-correcting, 
and so should not pose a major barrier in entering into the inhation phase. For 
F > H, the severity of this problem diminishes further since this damping term 
is an additional ehect that helps suppress the initial kinetic energy. Another sec¬ 
ondary detail is that above we only treated concave potentials. For double well 
potentials, V ~ (0^ —m^)^, as in the case of new inhation |p.l4|| , inhation requires 
the held to be well localized at the top of the potential hill 0 = 0 and almost 
at rest. These requirements are necessary, otherwise as found in various studies 
p.19 ], the held can easily break up into several small domains with randomly 
varying signs of the amplitude. For the supercooled case, new inhation, the basic 
conclusion about the initial conditions required for inhation has been that they 
are not very robust ||119|, |120| , |107|| . The inclusion of a dissipative term F > H 
will not necessarily improve this situation. On the one hand, such a term helps 
considerably to damp kinetic energy, thus allowing the inhaton to remain at the 
top of the hill longer and drive inhation. On the other hand, if such a damping 
term arises too early before inhation is to begin, since it freezes the evolution 
of excited inhaton modes, it could prevent the initial inhaton held conhguration 
from equilibriating to 0 = 0 . 

Finally in the above discussion F 7 ^ 0 generally represented warm inhation 
dynamics, although even for supercooled inhation, if radiation is present at the 
onset of inhation, such as in the new m and thermal || 121 || inhation pictures, 
the ehective evolution of the inhaton could have a damping term of the form 
F0. Since our focus is on the initial phase at the onset of inhation, it is possible 
the ehect of such damping terms also may be applicable to the initial condition 
problem in some supercooled inhation models. 


6.3.2 Radiation component 

So far we only considered the inhaton held system, but in addition the universe 
could possess some background component of radiation energy density, pj?®. To 
realize inhation, minimally it requires the vacuum energy density to dominate 
in the patch, pj?® < Pv If initially pj?® is larger than pv, expansion of the universe 
will red-shift the radiation. Thus provided the vacuum energy sustains itself 
during this early period, inhation eventually will start. This is the standard new 
and warm inhation pictures. 

One could imagine that in some small patch inside a larger causally created 
spatial region, the inhaton held conhguration is reasonably smooth and is sus- 




























taining a sizable vacuum energy of magnitude 


Pv 



( 6 . 22 ) 


As the universe expands pj?® decreases and provided pv sustains itself, inflation 
will begin in the patch. However at this moment, there could be large energy 
fluctuations in the radiation bath or other helds that are just outside the putative 
inflationary patch. One would like to understand how probable it is for such a 
situation to curtail the inflation that had started in the patch. 

To answer this question based solely on causality considerations, one should 
consider the case of a perturbation starting on the patch boundary and moving 
across the patch at the speed of light, and make the extreme assumption that 
as the perturbation overruns regions of the inflation patch, those regions convert 
back to being non-inflating. The question is what minimal initial inflationary 
patch size is needed so that its expansion under inflation is faster than its con¬ 
traction due to the impinging perturbations. For the case of flat geometry H = 1, 
this question was addressed in ||1U7|| and they concluded that the patch size should 
be at least 3 times the Hubble radius in order for inflation to succeed in enlarg¬ 
ing the patch. For a general H, we can address this question by evaluating the 
maximum distance a perturbation can travel in the patch between the initial 
and hnal tf time of inflation 


-q 


a^pert = a{t) / a ^dt 


(6.23) 


Note that since a is growing rapidly, there is negligible error in taking tf ^ oo 
which implies a^pert essentially is the event horizon. We have examined Eq. ( |6.23| ) 
for a variety of supercooled and warm inflation models and generally hnd Xpert ~ 

^mas* 

Thus in the most ideal case, an inflationary patch smaller than the event 
horizon can not be stable. Here an important point of syntax should be noted, 
that the generic size here is the event horizon and not the Hubble radius, although 
for the flat case, H = 1, both are of the same order. The above is the ideal bound 
based on causality. However, realistic dynamics also must be considered. 

One case which corresponds to external perturbations entering into an ini¬ 
tially inflating patch is the case treated in Subsect. HIA of the mixing of high 
wavenumber modes. For the inflaton held interacting with itself or with other 
helds, we argued above that if a dissipative term of the form F^ is present in 
the inhaton ehective evolution equation, then mode mixing up to wave-numbers 
k/a < r will not occur within the time scale relevant to the initial condition 
problem. 

In terms of the radiation bath, suppose a small vacuum dominated inhation¬ 
ary regime emerges, which is immersed inside a larger region containing radiation 
and gradient energy density. Since the inhationary patch has negative pressure. 


69 






the surrounding radiation will flow into it. The degree to which the pressure 
differences are signiflcant to this process depends on the magnitudes of the radia¬ 
tion, gradient and vacuum energies in both regions as well as detailed dynamical 
considerations [|lTl| , |122|| F|. Nevertheless, suppose the radiation tries to become 
uniform over the entire region, including the patch. Since the patch started to 
inflate, it meant that initially pv ~ Pr®- As time commences, provided pv re¬ 
mains constant, since expansion of the background will decrease it means 
the amount of radiation that flows into the inflationary patch will be less than 
Pv. Thus irrespective of whether the inflaton dynamics is supercooled or warm, 
the patch should sustain a type of warm inflation due to this influx of radiation 
from the regions that surround it. 

A thorough understanding of the initial condition dynamics in presence of a 
background radiation component is complicated. To our knowledge, no quantita¬ 
tive analysis has been done along these lines. The review ||107|| only illustrated the 
problem with the example quoted above Eq. (|6.23|) based on causality considera¬ 
tions but offered no dynamical examples with respect to a background radiation 
component. Here we have offered one scenario where the background radiation 
component should not prevent inflation from occurring. In particular, this exam¬ 
ple demonstrates that causality bounds on the rate at which radiation or other 
energy fluctuations enter the patch are not the only aspect of this problem. In 
addition, the magnitude of the entering radiation must be adequately large to 
overwhelm the vacuum energy. We have offered arguments above that indicate 
this latter requirement is not generic. 


6.4 Dynamical effects on embedding 


As discussed in Sec. |0| , a causally favourable embedding requires the inflationary 
patch to start smaller than the background MAS and then grow larger than 
it. A minimum condition for this to happen is that the particle horizon of the 
background eventually becomes greater than the background MAS. 

For p = up, the particle horizon is given by 


d{t) = 


dx 


H 


'o -|- ^{t)x 




(6.24) 


Comparing this equation with Eq. (|6.12|) it can be seen that the particle horizon 
will only become greater than the MAS if a; < 1/3. This condition was also noted 
in ||123|| from a slightly different perspective. There the patch was taken to start 
larger than the background MAS. 

^The dynamics involves viscosity effects at the interface between the inflation and back¬ 
ground regions. This viscosity is unrelated to the dissipative coefficient T in Sec. 3.3.1 , which 
represents the damping of the scalar field amplitude within the inflationary patch. 
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It follows that our proposed embedding will not work for a pure radiation 
background which has lo = 1/3. However, in general the background will consist 
of radiation and an inhomogeneous scalar held. As discussed in Sec. |6.3|, the 


gradient energy of the scalar held has a; = —1/3 and the potential energy has 
a; = — 1 while the kinetic energy has uj = 1. It follows that if the background 
is gradient dominated our causal embedding scheme will be viable, whereas our 
scheme fails if the kinetic term dominates. However, the kinetic energy of the 
scalar held will be suppressed if the dissipative coefficient V of Sec. is suffi¬ 
ciently large. This is already required inside the putative inhation patch in order 
to realize inhation. Thus, it is not unreasonable to expect T in the background 
region to be of the same magnitude. 


6.5 Conclusion 


This Chapter has investigated the initial condition problem of inhation from 
the perspective both of spacetime embedding and inhaton dynamics. Our study 
has highlighted two attributes of this problem which have not been addressed 
in other works. First, from the perspective of spacetime embedding, we have 
observed that the global geometry can play an important role in determining 
the size of the initial inhationary patch that is consistent with the weak energy 
condition. Second, from the perspective of inhaton dynamics, we have noted that 
a.T(f) damping term could alleviate several problems which traditionally have led 
to large scale homogeneity requirements before inhation. 

The purpose of this Chapter was to note for both these attributes, their salient 
features with respect to the initial condition problem. In the wake of this, several 
details emerge that must be understood. Below, we will review the main result we 
found for both attributes and then discuss the questions that must be addressed 
in future work. 

For a causally generated patch a successful embedding can be achieved if 
the patch does not contain an anti-trapped region. We have shown that the 
MAS size can be arbitrarily larger than the Hubble scale provided H is made 
small enough. So if the patch Hubble horizon is taken as the minimum stable 
patch size, then the patch does not have to contain an anti-trapped region if 
Hinf < 1. This generalizes the analysis of which was only for H = 1. However, 
without the effects of damping, it appears that the event horizon, not the Hubble 
horizon, is the minimum stable patch size. For the de Sitter case one can see 
analytically that the event horizon is equal to the MAS size regardless of H and 
numerical calculations suggest the same is true for power law inflation. However, 
radiation damping of perturbations could stabilize a patch smaller than the event 
horizon. In this case an open geometry for the patch would allow the patch not to 
contain an anti-trapped region and thus allow a causal embedding in an expanding 
background which does not violate the weak energy condition. Eventually the 
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patch should develop a MAS within it, but by then it could have expanded to be 
larger than the background MAS. 

For the dynamics problem, with respect to the scalar held the new consider¬ 
ation was the effect of a F^ damping term. We found that such a term could 
suppress many of the effects from initial inhomogeneities of the inffaton, which 
in studies traditionally done without this term lead to important impediments to 
entering the inflation regime. It appears evident that inclusion of such a damp¬ 
ing term will lead to qualitative differences in the initial condition problem. The 
most interesting outcome is initial inflationary patches smaller than the Hubble 
radius 1/H may be able to inflate. 

This Chapter examined the consequences of damping terms but did not delve 
into their fundamental origin. For the cosmological setting, such damping terms 
are typically associated with systems involving a scalar field interacting with fields 
of a radiation bath. In this case, such damping terms have been found in first 

EH 


principles calculations for certain warm inflation models 


although more 


work is needed along these lines. It is worth noting here that in the early stages 
of certain supercooled inffaton scenarios where radiation is present, in particular 
new [ 114 | and thermal [1^] inflation, a careful examination of the dynamics may 
reveal damping terms similar to this. Since the initial stages are the crucial 
period for the initial condition problem, if further study supports the importance 
of such damping terms, it may be useful to better understand damping effects 
also in such scenarios. 

Since the most suggestive situation for the damping terms is where in addition 
a radiation component is present in the universe, in Subsect. IIIB we also studied 
the effects of this component on the initial condition problem. Specifically we 
studied the case most suggestive for the initial stages of new and warm inflation, 
where a small inflation patch is submerged inside a larger radiation dominated 
spatial region. Our main observation has been that the minimal condition for 
the patch to inflate is that its vacuum energy density must be larger than the 
background radiation energy density. Provided the vacuum energy sustains itself, 
since expansion of the universe will dilute the radiation energy density, it is not 
clear-cut that the external radiation energy can act with sufficient magnitude to 
impede inflation inside the patch. 

We have shown that the gradient terms, due to the inhomogeneous scalar 
field in the background space time, make it possible for the patch boundary to 
overtake the background MAS. However the kinetic energy of the scalar field 
must not be dominant for this to happen. This can be ensured by including a 
damping term in the scalar field equation. 
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Chapter 7 
Conclusions 


In this Chapter the results of the thesis are summarised and current and possible 
extensions of the work are discussed. 

In Chapter adiabatic and entropy perturbations from multi-field models 
were studied |l[]. A general perturbation was decomposed into components paral¬ 
lel (adiabatic) and perpendicular (entropy) to the background trajectory. We de¬ 
rived the evolution equations for the entropy and adiabatic perturbations. From 
these equations it was evident that on large scales the evolution of the adiabatic 
perturbation is driven by the entropy perturbation, whereas the entropy is not 
driven by the adiabatic perturbation. We found that inflationary potentials that 
have a small effective entropy mass can generate non-negligible entropy perturba¬ 
tions. Also, it was shown that it is only for a non-straight background trajectory 
that the entropy will source the adiabatic component on large scales. 

We also analysed how correlations in the adiabatic and entropy components 
can arise in multi-held models. We found that when the adiabatic component is 
driven by the entropy component, it becomes correlated with it. The magnitude 
of this correlation depends on the relative magnitude of the homogeneous and 
inhomogeneous solutions of the adiabatic evolution equation. The magnitude of 
the isocurvature perturbation, in the radiation era, depends on the magnitude 
of the entropy perturbation at the end of inhation and the details of how the 
transition from the inhation to radiation era takes place. But, on large scales, it 
is independent of the adiabatic perturbation. 

We analysed the non-interacting double inhation model using our decompo¬ 
sition. Our approach was useful in identifying the reason for the features of the 
resulting power spectra during the radiation era. 

The question of whether preheating could ahect the spectrum of perturbations 
on scales relevant to structure formation was addressed in Chapter 0 0, §• A 
range of models that might have this property were examined. A potential which 
has efficient preheating and does not suppress the entropy perturbation during 
inhation was shown to be necessary for preheating to ahect large scale structure. 

We found that when the mass of the entropy held is larger than the Hubble 
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parameter during inflation, the entropy perturbation will be exponentially sup¬ 
pressed and so preheating will not be able to effect large scale structure. We 
also found that our new entropy evolution equation was much more numerically 
robust, in tracking this suppression, when compared to using the usual held per¬ 
turbation equations and constructing the entropy perturbation algebraically. 

The formalism discussed in Chapter ^ (which originally appeared in [jl[) has 
subsequently been extended by Bartolo et al. [^. They derived generic slow roll 
expressions for the adiabatic and entropy perturbation. One question that we be¬ 
lieve still needs to be resolved is how the entropy perturbation is transferred from 
the inflation to the radiation era. When there are couplings between the scalar 
fields, how these couplings influence the transfer is still not properly understood. 

In Chapter ^ we studied the effect of including correlated adiabatic and en¬ 
tropy perturbations on parameter estimation using the CMB [Q. We proposed 
a generic form for the initial power spectra which should encompass a broad 
range of two-field inflation models. We then modified the program CMBFAST 
and showed how considerable savings in computation can be achieved by using 
library spectra. 

We found that the current CMB data was consistent with a correlated CDM 
entropy perturbation of as high as twice the magnitude of the adiabatic perturba¬ 
tion. We also found a degeneracy between the primordial slope and the amount 
of correlated CDM entropy perturbation. 

Subsequently Bartolo, et al. ||124|| have used the formalism we developed in 
[Q to derive consistency relations between the scalar and tensor slopes for two 
held models. These consistency relations could provide a detectable signature for 
two-field inflation models in the CMB. 

In Chapter ^ we studied the the effect of having an extra spatial dimension on 
density perturbations in the brane world scenario 0. We analysed super-Hubble 
scalar perturbations using a covariant approach, and found that they evolve dif¬ 
ferently from general relativity due to bulk effects. The metric perturbation is no 
longer constant on super-Hubble scales during inflation and high energy radiation. 

We examined how the bulk Weyl curvature can be treated as an additional 
fluid on the brane with possibly non-adiabatic perturbations. We showed that 
on large scales the density perturbations on the brane are described by a closed 
system of equations, while on small scales there is no equation describing the 
evolution of the anisotropic pressure of the projected Weyl curvature on the 
brane and so the full five dimensional problem has to be studied. Subsequently 
Langlois et al. have shown how even the SW effect cannot be computed on 
the brane without solving the 5D equations. 

In Chapter we examined how non-linear inhomogeneities effected the start 
of inflation |p. We extended previous work to show that even for a non-flat 
geometry it was not possible to initiate inflation if there were perturbations on a 
scale smaller than the Hubble horizon. 

We then examined the effect of including a damping term due to decay of 
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the inflaton into radiation during inflation as in the ‘warm inflation’ model. We 
found this damping term set a new minimum wave length for the non-linear 
perturbations that could be smoothed by inflation. This minimum wave-length 
could be smaller than the Hubble distance and so this allows a causal mechanism 
to be responsible for the initial configuration of the inflationary patch. 
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